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Abstract: Towards the end of brane inflation, the brane pair annihilation produces 
massive closed strings. The transfer of this energy to Standard Model (SM) open string 
modes depends on where the SM branes and the brane annihilation are located: in the 
bulk, in the same throat or in different throats. We find that, in all cases as long as the 
brane annihilation and the SM branes are not both in the bulk, the transfer of energy 
to start the hot big bang epoch can be efficient enough to be compatible with big bang 
nucleosynthesis. The suppression of the abundance of the graviton and its Kaluza-Klein 
(KK) thermal relics follows from the warped geometry in flux compactification. This 
works out even in the scenarios where a long period of tunneling is expected. In the 
multi-throat scenario, we find a dynamical mechnism of selecting a long throat as the 
SM throat. We establish three new dark matter candidates: KK modes with specific 
angular momentum in the SM throat, those in the brane annihilation throat, and different 
matters generated by KK modes tunneled to other throats. Since the latter two couple 
to the visible matter sector only through graviton mediation, they behave as hidden 
dark matter. Hidden dark matter has novel implications on the dark matter coincidence 
problem and the high energy cosmic rays. 



Contents 



|T|. Introduction and summary 

|P| Introduction 
|1.2| Summary 



KK modes in GKP compactification 
|2.1| KK modes in a throat 



27| Tunneling between two throats 



273] Resonance effect in tunneling 
|3|. KK interactions 

|3j. Production of KK modes from brane annihilation 
|5|. Single throat heating 
|6|. Double throat heating 



6.1 KK evolution in A-throat 



6.2 Tunneling from A-throat to S-throat 



673| KK evolution in S-throat 
0. The bulk case 
§|. Multi-throat heating 



1] Throat deformation in Hubble expansion 
[8^ Selection of long throats 

|8.3| Warped KK dark matter and hidden dark matter 
Hidden dark matter 

Q| Hidden dark matter and dark matter coincidence problem 



972] Hidden dark matter and cosmic rays 
9T3| Comparison to other dark matters 



10. Conclusions and discussions 



[A], a KK cannot decay to all gravitons 



1. Introduction and summary 



1.1 Introduction 

One of the most interesting aspects of brane inflation [I], |3] is its close connection 
to experiments. This includes the spectral index, tensor mode, non-Gaussianities and 
cosmic strings — ranging from astronomical observations, gravitational wave detections 
to cosmic microwave background radiation (CMBR) measurements. 

Flux compactification in type IIB string theory has become an ideal setup to real- 
ize the brane inflation. This is not only because it provides promising mechanisms to 
stabilize moduli |J §, || and therefore a consistent string theory setup, but also because 
the flux induced warped space has become a very interesting new ingredient for model 
building. 

First, warped throats provide different mechanisms to achieve brane inflation. Anti- 
D3-branes naturally settle down in throats. The warped geometry red-shifts the brane 
tension and reduces (easily by a large factor) the mutual attraction between the branes 
and antibranes. This effect helps the slow-roll inflation, as in the KKLMMT scenario 
[0, || . Since the inflaton (the position of a D3-brane) is an open string mode, its kinetic 
term appears in the form of a DBI action. In warped space this imposes a speed-limit 
to slow down the brane motion in case the radiatively corrected potential becomes too 
steep. This mechanism generates the DBI inflation [§, [K| [IT], [12]]. Warped geometry 
together with the DBI action essentially guarantees enough e-folds of inflation. At the 
same time, effects of these two different mechanisms, and its combinations [13], can give 
rise to interesting observational signatures. 

Second, the simple geometric picture allows us to explore variations of the simplest 
scenario involving multiple throats and branes. The multi-throat configuration arising 
from the flux compactification provides a natural setup to generate candidate inflatons 
pT}| . Antibranes settle down in throats and annihilate against fluxes quantum mechani- 



cally [14], In 4-d view point, this process can start the brane inflationary epoch by 
eternally creating bubbles in the old inflation background. Within each bubble, the flux- 
antibrane annihilation naturally generate many D3-branes. The vacuum energy density 
within the bubble can still be large. These liberated D3-branes can roll out of the throat 
to other places and generate different types of brane inflation, such as the slow-roll or 
DBI inflation, whenever required conditions are satisfied. Different bubbles can also cor- 
respond to annihilations in different throats - the frequencies of which are determined 
by the relative lifetimes of the antibranes. In some cases, multiple branes may also be 
used to dynamically tune a flat potential |16 . 



Third, hierarchical range of string scales in multi-throat compactification provides 
different opportunities to observe stringy physics. Low string scale throat can make 
strings relatively easier to be observed in colliders ||17||, or can imprint stringy information 



on the density perturbations at large scale [12, 18]. High string scale throat can store 
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high tension cosmic strings at the end of brane inflation. The allowed range of cosmic 
string tension overlaps with observational abilities |19[ ^(| ^T], p2| . In this paper, we 
propose another possible consequence of brane inflation, namely hidden dark matter. 

The crucial step that links the inflationary epoch to the hot big bang epoch is 
the heating at the end of the inflation. This is known as the graceful exit. However, 
the heating in both the single- and multi-throat configurations remain so far the least 
understood aspect of brane inflation. Namely, how can inflationary energy be efficiently 
transfered to heat up the Standard Model particles, and be compatible with the well- 
understood late-time cosmological evolution? This is the heating problem (also called 
reheating or preheating problem). To see why this is quite a non-trivial question, we 
first look at the end process of the brane inflation. 

The brane inflation typically ends when branes annihilate antibranes. Significant 



insights have been gained into such a process recently [B3, 24, Go]. Tachyonic modes ap- 



pear when the brane- ant ibrane distance approaches the string scale and the annihilation 



process may be described by the tachyon rolling |26|, |27| . No matter whether there are 
adjacent extra branes surviving such an annihilation, the initial end product is expected 
to be dominated by non-relativistic heavy closed strings [28, 29, 30 1. These will further 
be involved to create lighter closed strings, KK modes, gravitons and open strings. (In 
this paper, we use KK modes as an abbreviation of the KK modes of massless closed 
strings.) We know from observations that, during big bang nucleosynthesis (BBN), the 
density of gravitons can be no more than a few percent of the total energy density of the 
universe. The rest is contributed by the Standard Model (SM) particles (mostly photons, 
neutrinos and electrons), which are open strings attached to a stack of SM ( ant i-) branes. 
We also know that the density of any non-relativistic relics can be no more than 10 
times that of the baryons. Therefore the question becomes how the brane annihilation 
products, originally dominated by the closed string degrees of freedom, can eventually 
become the required light open string degrees of freedom living on the SM branes, with 
a negligible graviton density and a non-lethal amount of stable relics. 

Depending on where the SM branes and brane annihilation are, we may consider 3 
possibilities : 

• when the SM branes are in the bulk (e.g., D7-branes wrapping a 4-cycle), while the 
brane annihilation happens in a throat; or vice versa; or they are both in the bulk; 

• the single throat scenario, in which case the SM (anti-)branes are in the same throat 
where the inflationary brane annihilation takes place; note that the model can still have 
other throats around; 

• the multi-throat scenario, where the inflationary branes annihilate in one throat (A- 
throat) while the SM branes are sitting in another throat (S-throat). 

Both the slow-roll and DBI inflation can be realized in different scenarios. Depending 
on whether one uses the throat to solve the hierarchy problem or to stabilize the cosmic 
strings, different models have different requirements. For DBI inflation, a single A-throat 
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may both fit the density perturbation and have a large hierarchy, because its prediction 
on the density perturbations can be independent of the A-throat warp factor []TT], |12| . In 
slow-roll brane inflationary models, the CMBR data typically requires only a modestly 
warped throat while solving the hierarchy problem requires substantially more warping, 
so the multi-throat scenario seems to be preferred as it can easily accommodate this 
very different warping properties. In both inflation models, the single throat versus the 
multi-throat scenario has another important implication. Brane annihilation is expected 
to produce all types of strings/defects that are not forbidden. This includes cosmic 
strings, which are simply massive large fundamental strings, D-l strings and/or axionic 
strings. They tend to survive much better in a throat without branes. In the single 
throat scenario, due to their interaction with SM branes, they are unlikely to survive 



cosmologically due to their potential instabilities [21, [31], B2L So the prediction 



of cosmic strings is much firmer in the multi-throat scenario, provided heating after 
inflation is not a problem. This calls for a careful examination of heating in the multi- 
throat scenario. 



A number of studies have been done to address this heating problem 
An important observation is that, because the KK mode wave function is peaked at 
the bottom of the throat, its interaction with particles located at the IR side is much 
enhanced compared to that with the graviton. This is essentially along the line of 



Randall-Sundrum (RS) model [17|. Because of this, the graviton emission branching 
ratio during the brane decay and KK evolution is suppressed by powers of warp factors 
^4, |3j| 36]. Such a suppression is absent if the brane annihilation happens in the bulk. 



But several concerns have also been raised |]S4" , SJ), SB]. Since the KK modes in a 
throat can have (approximately) conserved angular momenta, it is important to estimate 
its relic density. This was previously estimated to be extremely large |SJJ, and therefore 
any long-living KK relics seemed to be very dangerous to the BBN. Serious problems 
have also been raised in the multi-throat scenario where the question of heating become 
especially sharper ]34 , [3~6]1 . For example, in the multi-throat heating, the energy released 
from the brane annihilation has to be transfered to anther throat by tunneling, which is 
relatively a very long process. During this period, it seems that the KK modes in the A- 
throat is doing nothing but annihilate to gravitons. Requiring the initial annihilation rate 
to be smaller than the tunneling rate imposes very restrictive conditions on parameters 
of the model |36j . If generalized to the warped compactification where the tunneling rate 
is even much smaller, the analogous condition would rule out the multi-throat heating 
scenario. 

However, a more detailed thermodynamic evolution of the heating process, especially 
that of the KK particles, has not been carefully studied. Such a thermal history is 
important for obtaining various estimates and assessing the viability of the brane inflation 
heating. This paper is a step towards this direction. Other new ingredients are, instead 
of the Randall-Sundrum (RS) 1-d warped geometry, we consider a 6-d warped geometry 
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more relevant to the brane inflation in IIB theory (although the results can be easily 
evaluated and generalized to the RS case). This will cause interesting modifications in 
various estimates on tunneling rates and interaction rates. In a realistic compactification, 
throats are typically separated in the bulk, which tends to generate resonance effects in 
the tunneling from one throat to another. We expect the compactification volume to be 
dominated by the bulk, another important ingredient in the success of the graceful exit. 



We will find many qualitatively different results compared to those in Ref.[34, 35. 36 



For example, the KK modes will decouple typically after they become non-relativistic. 
Because of the warping enhanced KK self-interactions, their relic density turn out to be 
suppressed by powers of warp factors and the bulk size. In the multi-throat case, the 
KK relic abundance is further suppressed due to an extra matter dominated tunneling 
phase. Also for the multi-throat case, although the graviton production rate is initially 
much greater than the tunneling rate, the rate gets diluted by the spatial expansion. So 
the graviton production is significant only in the early epoch and its eventual abundance 
is also suppressed by similar factors. The most stringent constraint on the tunneling is 
just that it cannot be too long to over-cool the universe. 

So a final picture emerges is that a tower of KK particles gradually lowers as the 
universe expands, sometimes intervened by a long period of tunneling, and eventually 
transfers most of its entropy to the Standard Model particles at the bottom of the 
spectrum. The KK relic and graviton density are both reasonably suppressed by various 
factors due to the warped compactification. 

During the course of this study, we will find many other interesting phenomena 
arising from the multi-throat heating. First, there is a dynamical process that selects 
a long throat to be heated. This is because the dense spectrum in long throats makes 
the level matching of the energy eigenstates, a necessary condition for tunneling between 
throats, easier to satisfy. This may provide a dynamical explanation of the selection of 
the RS type warp space as our Standard Model throat in the early universe. Second, we 
find several new dark matter candidates. For example, we propose a novel type of dark 
matter - the hidden dark matter - which is generated during the tunneling heating. The 
hidden dark matter has many unusual properties compared to the usual dark matter 
candidates. 

Throughout the paper, we shall see the importance and richness of the warped space 
and multi-throat compactification appearing in many places, making them an even more 
attractive cosmological setup, in addition to the various aspects already mentioned at 
the beginning of this introduction. 

To simplify the analysis, we ignore closed string modes other than the gravity sector. 
The lightest modes in the other relevant closed string sectors (such as moduli) may 
either be more massive and stabilized throughout the whole process; or become part of 
massive particle spectrum with similar properties, and their inclusion will not make a 
qualitative difference on the basic issues. Similar situations apply to the gravitino and 
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its KK modes, so it is reasonable to expect that the gravitino density will not pose a 
serious problem cosmologically in this setup. In general, the qualitative properties of a 
simplified warp metric is very close to that of a more realistic metric, such as that of the 
Klebanov-Strassler (KS) throat f38fl , whose KK properties have been studied |39|, (40| . 
The simplified metric of S 5 allows a simple analytic treatment so much of the physics 
is more transparent. In instances when the qualitative properties of the metric of the 
warped deformed conifold differ from that of S 5 , we shall point that out explicitly. 



1.2 Summary 

To make reading easier, here we summarize by giving a brief history before the big bang 
nucleosynthesis according to the brane inflation heating in the multi-throat scenario. For 
simplicity, here we only keep the most important factors such as the warp factors hA,s, 
the length scales of the throats Ra,s, the string scale m s , the string couple g s and the 
Planck mass Mp\. The subscript A stands for the A-throat where brane annihilation 
towards the end of inflation takes place, while the S-throat is where the SM branes sit. 

Immediately after the inflation at t w gl^h^Mpi/m 2 ., all branes and antibranes 
annihilate in the A-throat and decay into heavy closed strings. The universe enters 
into a period of matter-dominated phase. These closed strings subsequently decay into 
gravitons and their KK modes. The branching ratio for graviton production is very small 
because the KK modes are exponentially peaked at the infrared tip of the A-throat, while 
the graviton, being the zero mode, is not localized. In fact, our estimate of the graviton 
production during this epoch is much smaller than that in Ref. |$5| due to the bulk size 
suppression. 

After a very short period At ~ g~ 2 h A 1 mj 1 , the decay ends and then quickly thermal- 
ize into a tower of relativistic KK particles at the bottom of the throat. The temperature 
is below the red-shifted string scale hATn s . The universe becomes radiation-dominated. 

As the universe expands, the temperature drops and the tower of KK particles one 
by one become non-relativistic and then quickly fall out of equilibrium. Because of 
warped space enhanced KK self-interaction, the abundance of these frozen-out species 
are suppressed by the warp factor and bulk size. 

This lasts until t ~ h~ A 2 R 2 A Mp\ when the temperature drops below the lowest KK 
mode, T ~ h A R~ A '■ The universe goes back to be matter-dominated. 

These lowest KK modes, which stable against decay to gravitons, enter a long period 
of tunneling process to other throats, in some cases going through the bulk resonance. 
During this period, the KK particles are decoupled from each other because of the red- 
shifting of the spatial expansion. The graviton production is negligible comparing to 
that in the early KK thermal history. Since the graviton density goes as a -4 while the 
KK modes go like a -3 (a is the scale factor), the graviton density rapidly becomes even 
more negligible. 
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Tunneling between two throats requires the energy level matching, and the chance 
of matching becomes more frequent for longer throats, since longer throats have denser 
spectra. Also because of this, the KK particles tunneled into these throats quickly decay 
into lower KK levels and no longer tunnel back. Due to these reasons, long throats 
are preferentially heated during the tunneling. It is reasonable to assume that the S- 
throat, at TeV scale, is most warped throat or among the most warped. So this process 
heats the S-throat, as well as generating possible new dark matter candidates in other 
throats. This long matter-dominated tunneling phase further suppresses the final KK 
relics abundance in the A-throat and the graviton abundance by delaying the Standard 
Model radiation-dominated phase. 

By t ~ ^tuni most energy has tunneled out of the A-throat. Here the tunneling rate 
T tun can be of order h\R~^ or h^R^ 1 depending on whether the bulk resonance happens 
or not. For example, in Ref. |7], |8|, [R|, matching brane inflation to the CMBR data, 
the warp factor Ha may range roughly from 1/4 to 10~ 3 . This will yield a reasonable 
tunneling time for BBN. In this case, the most heated throat turns out to be the S- 
throat. The KK particles in this throat thermalize to a much lower temperature and 
bring the universe back to be radiation-dominated. At this moment, the temperature in 
the S-throat can be above or below the local red-shifted string scale depending on the 
details. 

A similar story then applies to the KK tower in the S-throat. The difference is that, 
after the temperature drops below the lowest KK mode at t ~ hg 2 RgMp\, the universe 
remains radiation-dominated. This is because all the KK modes in the S-throat finally 
release most of their entropy to the Standard Model particles on the branes at the bottom 
of the mass spectrum. The KK relics is similarly suppressed as those in the A-throat by 
the bulk size and (much smaller) warp factors. This connects to the BBN. 

The KK modes in the bulk will tunnel to the throats, mostly because, without 
warping, they are very massive, and their wavefunctions spread out throughout the 
bulk. So we do not expect any KK relics in the bulk for the resonance case. We can also 
easily replace the S-throat with a bulk SM brane and the similar story will apply after 
KK tunnels out of the A-throat. 

Our calculations suggest that KK modes may supply three types of new dark matter 
candidates. 

First, since the wave functions of KK modes with specific angular momentum in 
the S-throat have no overlap (or very suppressed overlap) with SM branes sitting at the 
bottom of the throat, they do not decay into SM particles. As a result, they can be 
stable enough to become dark matter. Note that these are gravitational KK modes and 
that the Standard Model particles do not have KK excitations. So they are like the 



warped KK modes but different from the KK dark matter suggested in Ref. 

[[43| , [44] , [ll|. They interact with the SM particles via the exchange of (unstable) s-wave 
KK modes. In this sense, they are not too different from other dark matter candidates, 
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such as KK excitations of standard model particles or lightest supersymmetric particles, 
though they can have very different production and scattering cross-sections. 

Second, KK modes in another long throat away from the SM branes can also be dark 
matter candidates. If there are branes in such a throat, the KK energy can be converted 
to open string modes, much like what happens in the S-throat. Since these dark matters 
interact with SM particles only via gravitons, they will be hidden dark matter. They can 
be detected via gravitational lensing and rotation curves, but not in collider experiments 
or underground dark matter searches. These dark mater can have novel implications on 
the dark matter coincidence problem. 

Third, KK modes with specific angular momenta in the A-throat can also be quite 
stable. Their tunneling rate is highly suppressed compared to the s-wave KK modes 
that are responsible for tunneling in heating. This dark matter is formed as thermal 
relics much like the first type of the KK dark matter. But since they are decoupled 
from the SM particles except for graviton mediations, they also belong to hidden dark 
matter. Their masses are estimated to be very high, so their eventual tunneling to the S- 
throat and decay/interaction can produce very energetic events. Since these high energy 
events originate within the range of our galaxy but from outside the SM branes, they 
can generate cosmic rays that violate the GZK bound. 

This paper is organized as follows. 

In Sec. ^ we study the wave functions of the KK modes in GKP type of string 
compactification. We also study the tunneling properties of the KK particles between 
different throats. We estimate the tunneling rates in different cases where the bulk 
resonance is present or absent. 

In Sec. we estimate the cross sections of various possible interactions in the same 
setup. In both Sec. |2| and [| some quantitative differences between the GKP type and 
5-d RS type setup are pointed out. 

In Sec. H we outline the end process of the brane inflation - the brane-antibrane 
annihilation - and its various products. This leads to the initial condition of our main 
study of the KK thermal history. 

In Sec. ^ [| [7|and|8|, we study the thermal history of the KK particles in single throat, 
double throat, bulk and multi-throat cases, respectively. Estimates on the KK relic and 
graviton density are made. Interesting phenomena including the throat deformation in 
the Hubble background, the selection of long throats, the warped KK dark matter and 
the hidden dark matter are discussed. 

In Sec. ||, we study more properties of the hidden dark matter, focusing on the dark 
matter coincidence problem and possible implications on cosmic rays. 

Sec. [H] contains conclusions and discussions on future prospects. 
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2. KK modes in GKP compact ificat ion 



In this section, we study the wavefunction and dynamics of KK modes in throats. Previ- 
ous studies on brane inflation heating have been focused on the geometry of the orbifolded 
RS-type |42|, [34] . In this paper, we are mostly interested in the generalized GKP-type 
compactification ||, ||, [|6|, p8fl , namely, multiple throats with IR cutoffs whose UV sides 



are connected to a six- dimensional bulk. We will first study the wavefunction and spec- 
trum of a single throat case, then the tunneling of KK modes between different throats. 

2.1 KK modes in a throat 

In the GKP compactification, a KS throat is induced by three-form RR and NSNS fluxes 
around a conifold singularity. Except for the deformation around the IR cutoff of the 
throat, the gravitational and RR background generated by fluxes are similar to that by 
D3-branes with the same D3-charge. Therefore the analyses of KK modes dynamics is 
similar to that in the absorption cross section of black holes which have been studied 
extensively in the past JIT), [I8|, [|9], 0, |5T], |52], The difference here is that we will 



need to impose a different boundary condition, since we prepare our initial state inside 
a throat. 

We take the metric of the throat to be 

ds 2 = A- 1/2 (-dt 2 + dx 2 + h^dx^dx") + A 1/2 {dr 2 + r 2 dn 2 5 ) , (2.1) 
where we consider a simplified warp metric 

A(r) = 1 + ^ . (2.2) 

The r is the radial coordinate in six extra dimensions, and R is the characteristic length 
scale of the warped space. The 4D KK modes arise as metric perturbations h^, whose 
s-wave we decompose as 

= ^(a;)</>(r) . (2.3) 

The leading order equation of motion for the s-wave KK modes is a Laplacian equation 
in the background geometry, 

where m 2 = —p 2 is the mass of the KK mode with ip^ oc e ipx , and we have defined 
p = mr. 

To solve this differential equation, we redefine 

<j) = p~ 5/2 4) (2.5) 
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and get 



d 2 15 (mi?) 4 



dp 2 4p 2 



. 



;2.6) 



This differential equation is of Schroedinger type with zero energy, and the effect of the 
curved background geometry on KK modes is a potential barrier. We will study the 
inner and outer regions of the barrier separately and match them. 
The outer region is p 3> mR, where Eq. ( [2.6|) is approximately 



d 2 



15 



+ 1 



. 



dp 2 Ap 2 

The solution is given by the Bessel functions 

m = p il2 {H P )+iN 2 {p)). 



(2.7) 



(2- 



Here is the place where we have imposed the boundary condition, that is, we only 
consider outgoing wave in the large p region, because this is the KK mode leaking out 
of the throat. We can expand the solution (|2.8| ) in two different limits. For p ^> 1, 



(2.9) 



This describes an outgoing spherical wave. For mR <p< 1, 



— i- 



3/2 ^ 8 



(2.10) 



The first dots indicates the higher order imaginary terms that are related to the iiV^p) 
term in Eq. ( |2.8|) . The leading real term in Eq. ( [2.10| ) comes from the J%{p) term in 
Eq. @. 

The inner region is p <C mi?, where we can approximate Eq. (|2.6| ) as 



d 2 



dp 2 Ap 2 



15 (mR) 4 



P* 







(2.11) 



Changing variables, one can also bring this equation to Bessel's differential equation and 
get 



<Kp) 



p 



1/2 



mR 



AJ 2 



m 2 R 2 
P 



m 2 R 2 
P 



(2.12) 



We can again expand this solution in two different limits. For m R <C p mi?, 
7 *™ 3 R 3 „,m 7 R\ 4p 5 / 2 „,p l/2 , 



5p 



3/2 



P 



nm 5 R 5 



mR' 



(2.13) 
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Matching fl2l3|) and (|2~Tol) through p ~ mi?, we can determine the leading orders of the 
coefficients A and B in Eq. (|2.13 ) as 



A = -^ m - 3 R- 3 

7T 

B = -—m 5 R 5 . 
32 



(2.14) 



Note that there is also an imaginary part in B. We ignore this contribution since our 
eventual goal is to calculate the leading tunneling probability, in which such an imaginary 
term only leads to higher order corrections. For p <C m 2 R 2 , 




P 



2n m 2 R 2 



2 2 

(A - iB)e l( ^~^ ] + (A + iB)e 



2 2 
■ / m R 57r \ 

p 4 ) 



(2.15) 



This describes the outgoing and reflected incoming waves inside the effective potential 
barrier. To the leading order, they have the same amplitude. The small difference caused 
by the coefficient B leads to the tunneling probability 



A-iB 



A + iB 



7T s„s 

2S mR 



(2.16) 



We can perform a simple check on this tunneling probability by comparing fluxes of 
the outgoing waves between the inner and outer regions. To do this, we list the leading 
behavior of the wave function 0, defined in fl2.5|) , up to an overall normalization, 



2d2 . 



p < m z R 



m 2 R 2 < p < 1 



.i!^ m -BirV 3/2 



7T 

— P 

7T 



3/2 



COS 



m 2 R 2 



5ir 
T 



p > 1 



p -5/2 e i(p-tf) 



(2.17) 
(2.18) 
(2.19) 



For p <C m 2 R 2 we have equal amount of outgoing and incoming wave, the flux of the out- 
going component can be calculated by multiplying the flux density j = ^(0* ut V0 O ut — c.c.) 
with the area of warped unit five-sphere, 7r 3 A(r) 5//4 r 5 , and the warped unit longitudinal 
volume, A(r)~ l . Up to an overall normalization, we get 



Flux 



7T 

2i 



i>*r 5 



d_ 

Or 



— c.c. 



512 m~ 12 R 



(2.20) 



The region m 2 R 2 <C p 1 is the extension of the effective potential, and the KK 
particles are tunneling. For p > 1, we only have the tunneled outgoing KK particles. 
The corresponding flux is 



Flux 



2tt 2 



(2.21) 
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Comparing Eq. (|2.20|) and (|2.21| ), we get the tunneling probability ( [2.16|) . Here we have 
only considered the s-wave, the tunneling probability of the higher partial waves will be 
suppressed by factors of mR. 

The mass quantization of the KK modes comes from the IR cutoff of the throat at 
po = mrQ. Such a boundary condition will translate into a constraint on the wavefunction 
( [2.17 ) or its first derivative at po- The periodicity of the trigonometric function then 
leads to A(m 2 R 2 /p ) = it. So the KK mass is quantized in unit of 



Am = nhnR 



(2.22) 



where 



h = r /R 



(2.23) 



is the IR warp factor. The more realistic KS throat has a more complicated smooth 
shape at the IR end, but we expect the qualitative behavior of the mass quantization is 



similar [40|. 



2.2 Tunneling between two throats 

KK modes in the bulk can also be absorbed by a throat. Such a process is well analyzed 



previously when people study the black hole absorption cross section [g7|, [48], f|9], |50], 
51, [52], [53]. Here the boundary condition is to require only incoming wave in the inner 
region. We list the leading behavior of the s-wave in three different regions, 



p > 1 : 
m 2 R 2 < p < 1 : 
p < m 2 R 2 : 



.32^2 



57r\ 



M 

--m- 5 R- 5 , 

7T 



2 ir 2 p~ 3/ V ( 



(2.24) 
(2.25) 
(2.26) 



The tunneling probability turns out to be the same as Eq. ( |2.16| ). This leads to the 
absorption cross section for the s-wave [I 



n 



7T 



-m 6 R* 



(2.27) 



Now we imagine two throats - A and X. The A-throat is where brane and anti-brane 
annihilate and KK modes are generated. These KK modes can leak through the bulk and 
tunnel into another throat - X-throat, which is originally empty. We assume that the 
A and X-throat are separated by a distance D, and the bulk size is L. Let Jia = ta/Ra 
and hx = rx/Rx be the respective IR warp factors. The tunneling between the two 



throats requires the matching of the mass levels in two throats, n^/i^^ 1 



nxhxRx 1 , 
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within the energy widths. Here tia and nx are the KK levels in the A and X-throat, 
respectively. 

We will consider two interesting cases. In the first case, consider L > D > mT 1 . 
The size of the bulk are large enough, so that the KK modes with mass m may exist 
in the bulk. In this case, we assume such a mass quantization condition in the bulk is 
satisfied for m. The KK modes will go in two steps: first they tunnel out of the A-throat 
but remain in the bulk, propagating and bouncing around (in the extra dimensions), 
then gradually scatter and tunnel into the X-throat. For the first step, because of the 
tunneling probability ( 2.16 ) and the fact that the KK mode in A-throat has a bouncing 



period 21i a 1 Ra, the tunneling rate is 



7T 2 



A^buik ~ Yg^^^i^A 1 • (2.28) 

For the second step, the absorption probability is given by the cross section 

Pbulk^X ~ -jj, ' (2.29) 

and the bouncing period of KK modes in the extra dimensions is L. So the absorption 
rate is 

a 7r 4 o ,o R 8 Y 

1p K T'^MIi ' < 2 ' 30 ' 



bulk-*X 



The tunneling time for the KK particles to first go from the A-throat to the bulk and 
then from the bulk to the X-throat is the summation of these two tunneling time. So 
the final rate is the smaller of these two rates fl2.28j ) and (|2.30|) , 

Y A ^x ~ Min (T A ^ b uik, r w ^ x ) . (2.31) 



Because in ( |2.3U| ), L > n A h A Ra, we have T bu i k ^ x ^a^a-^x/^a- So the tunneling 



rate for ha ~ 1, Rx ~ Ra is roughly 

Ta^x < h\R~ A l . (2.32) 

The above analyses assumes that the energy of the KK modes in the A-throat satisfies 
the mass quantization condition of the bulk, namely the bulk KK mode also satisfies the 
mass matching condition as that for the two throats. This condition is the resonance 
effect that we will discuss in the next subsection. 

The second case is that such a bulk mass quantization condition is not satisfied. 
This also includes the case where the bulk size is too small, m < L~ x . Assuming the 
A and X-throats are still well separated, D ^> Ra,Rx, the KK particle has to tunnel 
two barriers at one time, since they cannot propagate in the bulk. Approximately, the 
tunneling probability will be the product of the probability ( |2.16|) for tunneling out and 
that for the absorption, and the rate for ha ~ 1 will be 

T A ^x ~ h\ 7 R A l . (2.33) 
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Of course one can always increase this tunneling rate by arranging special configuration 
of throats, for example bringing them closer and making the potential overlap. We will 
focus on the generic case. 

2.3 Resonance effect in tunneling 

There is an equivalent way of looking at the same problem involving the tunneling 
resonance effect in the presence of the bulk. Consider the same setup in the previous 
subsection. There are two barriers that the particle must tunnel through: the region 
between the A-throat and the bulk, r\ < r < and the region between the bulk and 
the X-throat, r 3 < r < r 4 . The KK modes have no barrier in the bulk, at r 2 < r < r 3 . 

The tunneling rate, or transmission coefficient, is straightforward to obtain in the 
WKB approximation. Beginning with the first barrier, the A throat to the bulk, the 
relation between the coefficients of the incoming waves and the outgoing waves is given 
by (see Merzbacher |54||) 



(2.34) 



6 A = 2exp^y dTy/V eff (T)j , (2.35) 

where T\ and r 2 are the classical turning point. Note that the effective energy in the 
analogous Schrodinger equation is zero. The variable r and the effective potential is 
defined in Ref. . The transmission coefficient from the A throat to the bulk is given 
by 

PA^buik = 4 [e A + — 

- A ■ ( 2 - 36 ) 

This WKB method amounts to solving the leading tunneling behavior of the differential 
equation ( p.5|) and (|2.6| ). Indeed, evaluating (|2.35| ), one obtains |L0|1 



P A ^buik ~ (n A h A f , (2.37) 

in agreement with the tunneling probability ( |2.16| ). 

Next we consider the probability of tunneling from the A-throat to the X-throat via 
the bulk. Up to some pre-factors, the matrix relating the coefficients of the incoming 
wave from the A-throat to the outgoing wave in the X-throat is given by 



X 

i(e A -ef) e A + e A l J\o e iW J \-i(e x -e^) e^ + e^ 1 



(2.38) 
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where W is the integral over the bulk 

W = r y/-V eff (T)dT . (2.39) 

For simplicity we assume Ra = Rx- Even though the throats have different warp factors, 
the WKB integral over the barrier between the bulk and the X-throat ends up to be the 
same as that from the A-throat to the bulk, Ox = ©a, since the KK quantization mode 

ua and nx obey n^/iA-R^ 1 = n xhxRx l ■ 

So the transmission coefficient from the A-throat to the X-throat via a bulk is given 

by 




(2.40) 



In the absence of the bulk, W = and T(A — > S) is very small, 

T A ^x ~ ©a 4 ~ (n A h A ) 16 . (2.41) 

This implies a tunneling rate 

T A -,x ~ n A 6 h A 7 R^ 

However, for KK modes satisfying 

W = (n w + 1/2)tt 

so cos W = 0, the transmission coefficient approaches unity 

T A ^x ~ 1 • (2.44) 

This is the well-known resonance effect. Notice that the location of the transmission 
peaks is determined by the same quantum condition for the bound states in the bulk. 
For very massive KK modes or large bulk, one finds that the resonances easily occur. 
This means that tunneling from the A throat to another throat passing through the 
bulk may not be suppressed at all if the state is just at the resonance point. However, in 
any realistic model, the initial KK wave packet will have some widths which we have to 
average over. Suppose h x ^ Ha- We expect a n^th KK mode will overlap with a small 
set of nxth KK modes so that the resonance effect will come into play. Or if some KK 
modes in the A-throat are in thermal equilibrium, then it is guaranteed that there are 
modes with the right energies to take advantage of this resonance effect. In the following 
we estimate the tunneling rate for these initial KK wave packets. 

For large G^, so that the penetration through the barriers is strongly suppressed, 
the transmission coefficient T has sharp narrow resonance peaks at these energies, with 
resonance width Ar. The distance between neighboring resonances is roughly 



(2.42) 



(2.43) 



Let us consider a wave packet that is localized in the A-throat at t = 0. To study 
the behavior of the transmitted wave packet near a resonance, we take the initial wave 
packet to have a mean energy E corresponding to a resonance. The packet has a width 
AE much bigger than the resonance width but much smaller than the distance between 
neighboring resonances, i.e., 

D > AE > Ar . (2.46) 

Around the resonance, 



cosW~ T [—) (E-Eo) (2.47) 



and sin W ~ 1 , so [ 54 



where 



For AE ^> AT, the total transmission probability for the incident wave packet is roughly 
AT / AE. One may interpret this result in the following way: the wave packet reaches the 
bulk in classical time. A fraction AT / AE of the packet is transmitted to the X-throat 
according to the usual exponential decay law with a mean lifetime 1/Ar. The rest of 
the wave packet is reflected. 

For the small band centered around the resonance point, using 

dW 

~ L > h~ A x R A , (2.50) 



dE 



we get the tunneling rate 



< n\h\R~ A l , (2.51) 



which agrees with ( |2.32| ). Note that this is effectively one barrier effect. This is because 
the resonance condition fl2.43| ) is just the KK quantization condition in the bulk, so if 
the resonance happens, the wave packet can stay and bounce around in the bulk. 

For a smooth distribution (e.g., a portion of a thermal distribution), the efficiency of 
transmission depends on Ar, AT / AE and the time it takes to replenish the transmitted 
part of the wave packet. Assuming the thermalization time is much shorter than the 
tunneling time, the rate of transmission is 



Ar 



.K^ y Ar<n>W, (2-52) 
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which agrees with that in Eq. (|2.42|) and ( |2.33|) . This can also be understood - a smooth 



distribution overall does not satisfy the resonance condition. 

In summary, for an initial KK wave packet in the A-throat (e.g. with — 1), if the 
bulk resonance happens, the particles can propagate in the bulk and the tunneling rate 
is enhanced for these particles, Ta^x ~ h A R^- Otherwise, without resonance or for a 
smooth distribution, the rate is much smaller, T^x ~ ^a-^-a ■ Notice that in either 
case, it is much smaller than T ~ h 5 R~ l for the RS 5-dimensional space-time case. 

3. KK interactions 

In this section, we study various cross sections of interactions among KK modes, gravi- 
tons and brane fields. Since we eventually only need an order-of-magnitude estimate of 
the thermal history, we will neglect the numerical factors for simplicity. 

The interactions among the KK particles and gravitons are determined by the 10-dim 
action 



Mf 



J d 10 Xy/=g^R w D Mf J d 5 tt dr A r 5 J d 4 xy/=glR 4 . (3.1) 



We will first express various couplings in terms of the 10-d Planck mass M\q and 
then convert it to 4-d Planck mass M-p\- From the action we can see that each 

throat contributes Mf R® to the 4-d Planck-mass-squared, while the bulk contributes 
Mf L e . Since the typical situation in the multi-throat compactification is that the size 
of the bulk L is bigger than the sum of length scales of all throats, ^ -Rj, we approximate 

Mpi « M± L 3 . (3.2) 

Expanding (|3.1|) in terms of the 4-d metric fluctuations h^ u , we get 

Mf J d 5 ttdr Ar 5 J d A x (dh ■ dh + h ■ dh ■ dh + h ■ h ■ dh ■ dh + ■ ■ ■ ) , (3.3) 

where we only write the relevant expression schematically by neglecting indices that 
should be properly contracted with rj^ and ignoring various numerical factors. After 
decomposing h^ v = ^2ip n ^u(x)(f) n (r) as in (|2.3|) and integrating out the wave functions 
fl2.17|) - ([2.19|) in the extra dimensions, we can get couplings for the 4-d KK modes ip^ v . 

To make the kinetic term for the KK modes canonical, we first integrate the extra 
dimensions for the quadratic term. It turns out that the dominant contribution comes 
from the tip region r ~ r , 

M* m~ 13 R- 7 hv 1 / d 4 x(di>) 2 . (3.4) 



where is the warp factor (|2.23|) . For the zero-mode graviton, the wave function in 



the extra dimension is a constant, for example, we will take it to be one. The quadratic 
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term is dominated by the bulk 



Mf L b J d 4 x(<9</> ) • (3.5) 

Therefore the canonically normalized fields are 

iP = M* m- 13 / 2 iT 7 / V /2 ^ > 

= M* L 3 ^o . (3.6) 

We next look at the interactions. We consider the case where all KK particles 
and gravitons have similar energy. This applies to the case of the thermal equilibrium 
that we will focus on later. So in terms of orders of magnitude, the 3-point interaction 
A3 J d 4 x ip ■ dip ■ dip corresponds to a cross section a ~ A|, and the 4-point interaction 
A4 J d A x ip ■ ip ■ dip ■ dip corresponds to a cross section o « X^E 2 , where E is the energy 
of the particles. We first look at the cubic interactions. 

Three KK modes: This is the KK self-interactions. Integrating along the throat in 
( p.3| ), the dominant contribution comes from the tip of the throat, where the KK modes 
wave function are peaked. In terms of the canonically normalized fields, we have 

M^R~ 3 ho 1 J d 4 x ip ■ dip ■ dip . (3.7) 

This corresponds to an interaction cross section 

(3-8) 



Two KK modes and one graviton: This is the KK modes annihilation process. The 
coupling of this interaction is still dominated by the tip region r ~ r , 

Mf 4 L- 3 J d A x ip ■ dip ■ dip , (3.9) 

where one of ip is ipo. This leads to an annihilation cross section 

* « ^ ■ (3-io) 

One KK mode and two gravitons: Because of the KK momentum conservation, 
namely 

dr r 5 A <f>(r) = , (3.11) 

such a cross section vanishes. This equation is in fact the orthogonality condition between 
the KK and the zero mode, so this conclusion can be easily generalized to more arbitrary 
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geometry or multi-gravitons. Since any compactification has a zero mode graviton with 
constant wavefunction ||0| , n-point function involving one KK mode and n — 1 gravitons 



is proportional to the orthogonality condition between the KK and graviton mode. Hence 
the KK mode cannot decay to all gravitons directly. This conclusion still holds when all 
other terms in the action are considered. The arguments are given in Appendix 
For higher KK mode, it can still decay to several lower KK modes plus graviton with a 
rate <C m 3 /Mpj. But for the lowest KK mode, the only way it can emit graviton is to 
annihilate or scatter another particle. 

Three gravitons: The graviton self-interactions are dominated by the bulk, their 
coupling and cross section turns out to be the same as the KK annihilation case, o ph 
1/M&. 

If branes exist at the IR tip of the throat, the open string fields on the brane will 
interact with KK particles and gravitons. These interactions are given by the following 
term in the action, 

J d 4 xd 5 ttdr A 1 ' 2 ^ ~ G" u d^Hd v H S(r - r ) <5 5 (ft - Q ) , (3.12) 

where G^ u = A 1 / 2 g liu , r^ and Oo are the radial and angular location of the D3-branes. 
The H denotes the brane fields, which for example can be a Standard Model Higgs. 
Expanding (ft.!2|) in terms of the 4-d fluctuations h^ u , one gets 



2r 2 



[ d 4 x (ifdpHdyH + hTdpHduH + h ■ h ■ OH ■ OH + ■ ■ ■) . (3.13) 

J V / r=r 



We denote the canonically normalized SM particle as 

H = —H . (3.14) 

The presence of the D3-branes breaks the conservation of the KK momentum. So for 
example, a KK particle can hit the brane and break into two SM particles. 

KK-SM interaction: Such a three-point interaction is determined by the second term 
in O, 

j ' d A x^ v d^Hd u H . (3.15) 



The cross section is 



'-i'BS- (3 - l6) 



If the KK particle is replaced by a graviton, we get 

1 



° « T72 ■ (3-17) 
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The quartic interactions can be analyzed in a similar way. In the following, we 
summarize both the cubic and quartic interactions according to the interaction type. In 
all cases, the cubic interactions are the most important. 



KK self-interactions: This includes KK + KK -> KK, KK + KK -> KK + KK, 
KK + KK — > + g and higher points interactions. The cross section of these 
three interactions are given by the three terms in 



pi 



respectively. Here n E parameterizes the particle energy scale, E = n^/io-R -1 - For 
relativistic species, n E > n KK , and for non-relativistic species, n E = ukk- In the 
last two terms, the relations R A ~ Ng s /m\ and Mpj ~ m^L 6 /g 2 have been used. 
These two terms are negligible to the first, so we have 



L\ 6 h7 2 



°KK ~ o 77^ • ( 3 - 19 ) 



Comparing to the RS case, we have the same warp factor enhancement 2 . In 
addition, the cross section in the GKP case is further enhanced by the bulk size 
L 6 /R e . 

Hard graviton production from KK annihilation: This includes KK + KK — > g, 
KK+KK — > KK+g, KK+KK — > g+g and higher points interactions. The cross 
section for the first reaction is (|3.10|) , and the third has an additional suppression 
factor E 2 jMl v The second process does not entirely annihilate the KK modes and 
was included in the previous case as a KK self-interaction. But this reaction emits 
a hard graviton and can be interesting as we consider the graviton production. The 

2 

cross section is tS-tt^-- So the first reaction dominates, 



N 2 M pl 

a K9 * ± ■ (3.20) 

Graviton interactions: This includes g-\-g—>-g,g-\-g—>-g + g,g + g^H-\-H, 
g + g — » KK + KK and higher points interactions. The largest cross section is 
that of the first, 

^ « 4r ■ ( 3 - 21 ) 

The latter three are suppressed by an additional factor of E 2 /M^ v So as usual the 
graviton interaction is Planck-mass suppressed. 
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KK-SM interactions: This includes H + H -> Kif, H + H KK + and 
higher points interactions. Their cross section is 





6 h- 2 


a 













(3.22) 



.f/ani graviton production from SM particles annihilation: This includes H+H — > g, 
H + H —>■ g + KK, H + H ^ g + g and higher points interactions. Their cross 
sections are 

n| E 2 



M 2 , 



(3.23) 



We note that the wavefunction Q2.17| )-( p.l9| ) used here is for s-wave. For higher 
partial wave, we will have a phase shift in the wave function and an extra angular part. 
The qualitative estimates we made in this section also apply to those KK modes with 
non- vanishing wavefunction at the bottom of the throat where the SM branes are sitting, 
such as those associated with S 3 in a realistic deformed warped conifold of the KS throat. 
There, the KK modes with S 2 angular momentum has vanishing wavefunctions at the 
origin, so they will interact with the SM particles via the exchange of the unstable s-wave 
KK modes. 



4. Production of KK modes from brane annihilation 

To study the heating process in brane inflation, we first look at the initial state from the 
brane-antibrane annihilation. We choose the time coordinate by choosing p ph M^/t 2 . 
So at the initial brane annihilation p ph h\T 3 , we have t « g l J 2 hJ^M-p\/m 2 . The 
string dynamics in brane annihilation can be described by a boundary conformal field 
theory [|23], ^4], |25|. Calculations in the tree level and loop level [f29], |3(J string 



amplitudes indicate that the dominant product is non-relativistic massive closed strings, 
no matter whether there are any branes left over from such an annihilation. These heavy 
closed strings will further decay to lighter states, such as lighter closed and open strings, 
gravitons and KK modes. The time scale for such annihilation and decay processes is 
of order red-shifted string scale, At ~ g^ 2 h~^m~ l , up to a constant factor depending 



on string levels and configurations [p5[ . The At is much smaller than to, so these decay 
happens right after t . 

As we will check later (in Eq. ( |5.5| )), it is easy to satisfy the condition that the spatial 
expansion of the universe is slow enough for these products to thermalize. So both the 
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decay and thermalization time is smaller than to. Therefore the initial energy density of 
the thermalized gas is approximately the same as the initial brane tension h\T 3 . For a 
gas of relativistic particles with g degrees of freedom, the thermalized temperature will be 
T ~ (gg s )~ 1 ^ 4: hA'm s . For g > g^ 1 this is below the red-shifted string scale in the throat. 
Hence these thermalized particles should be the graviton KK modes if T > HaR^^ as 
well as the brane fields (zeroth open strings modes) if there are branes left. 

Gravitons are also present as a decay product. Since their cross section is of order 
Mpj 2 , they are never in thermal equilibrium as long as time is larger than the Planck 
scale. So it is already important to estimate its initial abundance. 

If the brane-antibrane annihilation are in the bulk, the graviton and the bulk KK 
modes couplings to the heavy closed strings are of the same order of magnitudes. Then 
the initial graviton density will be a sizable fraction of the total energy density. Since 
during the BBN, graviton density cannot exceed a few percent of the total energy density 
of the universe, a more detailed analyses is required to determine this fraction. The 



analyses in Ref. [36] indicates that this is problematic. (However see a solution at the 
end of Sec. [7[) 

The situation is very different in the warped space. This is because, unlike the zeroth 
mode graviton, KK modes are peaked at the tip of the warped space. In the 5-dim RS 
setup, their coupling to the heavy strings located at the IR end of the warped space will 
be enhanced by a factor of h^ 1 relative to the graviton. Hence the graviton production 
is suppressed by a factor of h\ |34], [35], |36|] . The case that we consider here in the GKP 
setup is essentially the same except for the quantitative details. From the analyses in 
Eqs. (|515D - (gT?D , we know that the suppression for the graviton production is greater 



by a factor of 

~)V (4.1) 

To get an order of magnitude feeling, we note that, in fitting CMBR data, R/L ~ 1/3 
to 1/100 and Ha ~ 1/4 to 1/1000 are reasonable for inflation in A-throat [[?], || If 
inflation takes place as the D3-brane is coming out of another throat, as in the IR DBI 
case PH , p~2[] , then the value of Iia is more relaxed. 

So in the presence of any reasonable warp factor, the dominant products after brane- 
antibrane annihilation, string decay and particle thermalization is a tower of relativistic 
KK modes with possible brane fields at the bottom of the mass spectrum. 1 The thermal 
history of such a gas will be the main focus of our next few sections. The study of the 
evolution of such states will eventually lead us to conditions at the time of BBN, and 
hence provides a test whether this is a viable heating. 



1 The initial KK can be absent if To < IiaR^ , so we only have relativistic brane fields. This simplified 
case only happens for very small N A < g in the single throat case. 
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5. Single throat heating 



In this section, we study the time-dependent thermodynamics of KK modes in single 
throat case. In this case, the D3-anti-D3-brane annihilation, and therefore their eventual 
product KK particles, are in the same throat as the Standard Model branes. Our goal 
is to give an estimate of the relic KK mode and graviton production from the thermal 
history. To keep it more general, in this and next section, we will first obtain results 
in terms of general cross sections and tunneling rate, and then evaluate them using the 
results in GKP type compactification in Sec. ^| & |3|. It can also be evaluated in the 
Randall-Sundrum setup. 

Initially the Universe is dominated by a gas of relativistic KK modes that are created 
from brane and string decay. We know that the evolution of the energy density during 
such a period follows 



If the KK particle self-interacting rate T(t) is bigger than the Hubble expansion rate 
H(t), the KK particles are in thermal equilibrium. We will check shortly that this is the 
case for the period that we are interested in. The density of the thermalized relativistic 
particles is determined by the temperature, 



where g is the number of degrees of freedom of all relativistic particles in equilibrium. 
These particles include the KK and SM particles and their possible super-partners. For 
simplicity, we will take g to be a constant of the same order of magnitude as the SM 
particles ~ 100. Variation of g may come if the KK tower is very high at the early epoch, 
or in the absence of the SM particles. We will encounter such cases in the next section. 
The difference caused by the variation in g is not essential in this section. 

The time dependence of temperature T then follows from Eq. ( |5.1| ) and ( p.2| ) . 



As the universe expands, the temperature decreases. Each KK mode consecutively 
becomes non-relativistic as T drops below its mass. 

At T nr w m = nKxhoR~ l , the temperature crosses the mass of the n^xth KK 
modes (including all possible degenerate states) and these KK particles are becoming 
non-relativistic. From Q5.3D , the corresponding time is 



H 30 y ' 



(5.2) 




(5.3) 



t 



nr 



9 



/2 M Pl m- 2 . 



(5.4) 
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The condition that we can use the thermodynamics for this KK mode until this moment 
is 



T/H « « 9 T'a KK t nr > 1 , (5.5) 

where n is the number density of this KK mode and the velocity v of the KK particle is 
approximated as 1. This is easy to satisfy. For example, for o KK in ( |3.i9| ), 



T/H « g^UKKR^M- 1 (^) h, 1 



(5.6) 



If g 1 ' 2 « 10, n** w 1, m s w 10 10 GeV, A^ /4 » 10 (i^ « g s N A mJ A ), and M P1 w 
10 18 GeV, we need (|) Z^ 1 > 10 8 . This can be easily satisfied e.g. by L > lOi? and 
^0 < 0.01. 

The density of these non-relativistic KK particle will drop exponentially with tem- 
perature 

n nr = g(^P) 3/2 e- m/T , (5.7) 



2tt 



and their entropy is released to the other relativistic species. The self-interacting rate of 
such a particle also drops accordingly and this species quickly becomes decoupled. After 
t nr , the decrease of the interaction rate T = n nr aKKV is dominated by the exponential 
decrease of n nr . So when the KK particle is decoupled, T « H(td ec ), we can estimate 
the leading order of the particle number density as 



KK 

''dec 



#(W)oV> -1 (W) ~ H(t nr )a^ K « j^Mp/mV^ . (5. 



Here the changes in H and v from t nr to t& ec are neglected for simplicity. 2 These changes 
will give further suppression factor which is logarithmically weaker. 

The whole tower of KK modes is gradually lowered in this fashion, until the temper- 
ature drops below the mass of the last KK mode h^Rr 1 . At this point, all KK particles 
are non-relativistic and release most of their entropy to the relativistic SM particles. 

If any KK modes are stable because they carry some conserved quantity, they become 
relics after decoupling. It is important to estimate their relative abundance to the other 
relativistic SM particles. Since non-relativistic KK particles red-shift differently from 
relativistic SM particles, this relative abundance will be amplified in later evolution. 



Using (5.7), we can estimate T{td ec ) as m {h\^—^j . So, the velocity v at td ec is v « y/T/r 



( 3 \ — -L/Z 
m nT - ) ' wmcn changes logarithmically slower than tlkk, consistent with our assumption. The 

/ 3 \ —3/2 

same happens to H . An overall correction is a factor of (In -^wwj , which is smaller than one and 



makes the in (5.8) smaller 
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This amplification lasts until the transition epoch from radiation domination to matter 
domination (RDMD), at which time the density of these non-relativistic relics cannot be 
more than that of the radiations, even if they account for all the dark matter. 

We now estimate this relative ratio. At td ec , the ratio of the energy density of the 
decoupled KK particle to the radiation energy density is 

KK \ 

P \ n dec m . 

(5.9) 



p tot J , gm 4 g x l 2 mM?\o KK 

Summing over all possible stable KK modes and taking into account of the red-shift, at 
the time of the radiation-matter transition, we have 

/ KK 

c,kk _ I P 

iL RDMD - I tot 
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p KK \ a(t 



RDMD j 



dec ) 



n KK \ / n tot (j. \ \ 1/4 

P \ IP Kpdec) 



P tot J tde , V Prdmd 



^ g-^M-lp~^ MD a- K \ , (5.10) 



riKK 



where Prdmd is the energy density at the time of RDMD transition. Using the cross 
section okk in (|3.19|) , we have 3 

^MD^T.^ijX^hl. (5.11) 

n KK ^ ' Prdmd 

For the single throat case, this ratio can be made very small if h is small. Since 
Prdmd ~ 1 e V> M P \ « 10 18 GeV, we need ^ 10 -27 . For example, this can be 

satisfied if L w 100i? and ho w 10 -8 . A few comments are in order here. First, we 
can also evaluate (|5.10| ) in the 5-d RS case. It is still suppressed by the warp factor 



h,Q but without the bulk size suppression. So the condition becomes more restrictive, 
h < 1CT 14 . Second, as we will see in the next section, ^Ir^ md will receive another 
suppression factor in the double throat case due to an extra matter-dominated phase. 
Third, we note that we have been interested in the relic density resulted from the thermal 
history. Whether these relics can be stable enough is also important and case dependent 
[ j35[ , p~5| . In this paper, we show that even the most restrictive case, namely the stable 
relics, is compatible with a successful heating. 



3 The cross section of the KK self-interaction can be enhanced by the KK-SM particles interactions 
( 3.22p . Especially when we have many different SM particles. Such an enhancement can further decrease 



the ratio. 
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We now estimate the graviton production. They are produced through the KK and 
SM particles annihilation. Once produced, they are always decoupled and red-shifted as 
radiation. Therefore their density is 

p9 rav {t) „ /" T(ti) n 2 {h) ^ ^ (512) 

where we estimate the graviton energy as T, which is the energy scale of the reactions 
involved. The n « gT 3 is the number density of the KK or SM particles. From ( |3.20|) 
and (|3.23|) we know that an g ~ a Kg- The initial time to is set by initial energy density 



released from the brane annihilation, Hq m m ^pr, as in Sec. |[ So we get 

P 9rav (t) « g; 1/A g 1/4 Ml ims a Kg h^ . (5.13) 

In this integration, we can see that the dominant contribution comes from the initial 
moment to over a period of order t$. After that, p 9rav just red-shifts as 1/t 2 . Therefore 
the faction of gravitons in total radiation is 

„grav ' 
4ot 



^BBN or RDMD =(^) « ^1/4^^ . (5 . U) 



Z 5 / BBN or RDMD 

Using the cross section ( |3.20| ), we have 

BBN or RDMD ~ 9s 9 ^0 , 1^0. 10J 

which is naturally very small. 4 We note that the ratio Q5.14Q is proportional to the 
graviton production cross section a Kg, but is independent of the KK self-interactions 
cross section okk- This is the difference between the thermal dynamics and simple 
reactions. That is, KK and SM particles are in thermal equilibrium and their density 
is determined by the temperature rather than the cross section, while the gravitons are 
never in thermal equilibrium. The total graviton abundance is the summation of ( |4.l| ) 
and ( ETT5D . 



6. Double throat heating 

In this section, we consider the heating process in double throat case. In this case, 
the brane-anti-brane annihilation happens in the A-throat and there is no brane left 
in the end. The Standard Model lives on the branes in the S-throat. So successful 
heating requires efficient energy transfer from KK modes in the A-throat to branes in 
the S-throat. We study this process in three steps: the evolution of KK modes in the 
A-throat, the tunneling between A and S-throat and the evolution of KK modes in the 
S-throat. 

4 Higher KK particle can also decay to several lower KK particles plus graviton with a rate <§; m?/Mpy. 
Within a time interval t nr w g^ 1 / 2 M-p\m~ 2 , the decayed fraction is <C g s 1 ^ i g^ 1 / 2 N^ 1 l i ^^h, much 



smaller than (5.15) 
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6.1 KK evolution in A-throat 

The thermal history of the KK modes in the A-throat is very similar to that in the 
single throat case. The main difference is that at the bottom of the spectrum there is no 
longer SM particles. So even if the KK tower starts not very high, the particle degrees of 
freedom will not remain constant. We write the relation between g and the temperature 
T in the following ansatz 

9 -»"(w) TwaI+i (i^) T " (fU) 

where the value of 7 depends on the degeneracy of the states. For example, for 5-d RS 
case, 7 = 1; for AdS^ x S 5 , 7 ~ 5; for AdS$ x T 11 , 7 « 4 to 5. Sometimes, as in the last 
step, we will also write R in terms of brane tension T3, where a ~ iVT+ 4 . Therefore the 
thermal density of the KK modes is 

ijvy+4 

pit) » 7— , (6.2) 

and the time dependence of the temperature follows 

r W «(Mtf)*p£)*. (6.3) 

The KK modes at each level ukk become non-relativistic when the temperature drops 
below m ~ ukk^aRa a ^ 

t Anr w ) 2 M P1 m-— w n^ x 2 R 2 A M Fl h~ A 2 

~ ^Adec • (6.4) 

The decoupling time tyirf ec for these KK modes follows shortly after t^nr- The main 
difference between this case and the previous single throat case is that, as the temperature 
drops below the last KK mode in the A-throat, the universe enters a period of matter 
dominated (MD) phase. The lowest KK mode stays in the A-throat waiting to tunnel. 
The tunneling time t tun ~ ^tun ls estimated in Sec. 0. 

The higher partial waves of KK modes take much longer time to tunnel, and if stable 
enough, can become non-relativistic relics. This relics density can be estimated at the 
time when they become decoupled. The corresponding number density is 

UAdec ~ — I « tAnr* KK ■ ( 6 - 5 ) 

&KK 

The relative density of these KK relics at decoupling is 

P KK \ n Adeem 

) W M 2 /t 2 • (6 ' 6) 

P S t Adec M P\/ T Adee 
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This ratio will be amplified during any radiation domination epoch that follows. As 
we will see, this will start from when the tunneling is done (~ t tun ) and lasts until the 
RDMD transition in the S-throat. So at RDMD, the relics abundance is 



/ n KK\ f n KK\ i/l/2 r l/2 

n KK / P \ ~ P \ M V\ Hun 

* RDMD ~ I ~~~JoT J ~ \ ~7~toT I 1/4 

J RDMD nt ,„\P y t AH „ r Prmun 



TlKK z Adec 



Prdmd 



riKK 



n KK M V\PRDMD°KK f 7^1 ■ ( 6 - 7 ) 



^tun 



The main difference between Eq. ( |6.7|) and ( 5.10|) is the last factor in Eq. (|6.7|) . This 



further suppression is due to the matter dominant phase from tAnr to t tun . In the single 
throat case, after the last KK mode becomes non-relativistic at t nr , the universe remains 
radiation dominated. In the double throat case, radiation dominance is delayed until 

ttun- 

The graviton production during this period can be estimated using Eq. ( |5.12| ). Sim- 



ilarly the integration is mainly contributed from the early epoch around to, where to is 

1 

4& 



determined by the initial condition Hq ~ ^ ~ oIjt-- We get 



pT v ~ 9 ; 1/4 a 1/4 h A m s M^K 9 ^ ■ (6.8) 
The relative graviton density is 

(grav \ / , \ 2/3 

- g 7 ^aV±m s M Pl v Kg h A . (6.9) 

PA J BBN or RDMD V l tun J 

The last dilution factor is again due to the additional matter dominated phase from t Anr 
to t tun . Using (|3.20| ) we see that this is naturally a very small ratio. 



6.2 Tunneling from A-throat to S-throat 

Because of the very small tunneling rate (|2.31|) and fl2.32j ), the tunneling of the KK 
modes in A-throat effectively happens well after the temperature drops below the mass 
of the lowest KK mode. Therefore the tunneling is dominated by these non-relativistic 
lowest KK particles. The only way for the lowest KK mode to emit graviton is through 
interaction with other particles. But during the tunneling, they actually have long been 
decoupled from each other and their kinetic energy has been largely red-shifted. (For this 
gas, we know that one can assign an effective temperature T oc a~ 2 oc t _4//3 if we keep 
a correspondingly changing chemical potential.) The energy density evolves similarly 
as in Eq. ( |5.1|) except for a different numerical factor 4/3. To get the total graviton 
production, we integrate all particles as in Eq. ( |5.12|) . Although the tunneling time 
scale tends to be very long as we see from Sec. |2|, the graviton production from the KK 



annihilation in this period is negligible comparing to (5J3), as we can see from a similar 
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integration as in Eq. (|5.12|) . This is because the spatial expansion keeps on red-shifting 
the KK density and velocity. 

So the only constraint on the time scale of this period is that the spatial expansion 
during this long period does not cool the universe too much. This requires that the 
energy density at the end of the tunneling, PKK{Uun) ~ ^v^luni be larger than the 
energy density during BBN, (1 MeV) 4 . This gives a lower bound on h A - For example, 
for the resonance case or the non- resonance case, 

Ttun ~ h\R? or T tun ~ h A 7 R~ A l , (6.10) 

we need 

h A > (R A ■ MeV 2 /M P1 ) 1/9 or h A > (R A ■ MeV 2 /M P1 ) 1/17 , (6.11) 
respectively. If for example m s ~ 10 10 GeV and N A 4 ~ 10, we need 

h A >10~ 4 or h A >!0~\ (6.12) 
respectively. If we want to incorporate TeV baryogenesis in the S-throat, then we need 

h A > 10- 2 or h A > 10- 1 , (6.13) 

respectively. These are reasonable for an A-throat. Note that for the 5-dim RS setup, 
the bound is much looser. 

The KK modes tunneled into the S-throat will become thermalized in a rate 

/ 1 \ — 2(7+3) 

lW*tun) « n(t tun )a KK « {hsRs 1 ) (M Pl T tun ) ^ a KK , (6.14) 

where the time dependence of the number density of the KK modes n(t) is obtained from 
formula same as Eq. flK.lfl and (573). As long as this rate is larger than the tunneling rate, 



which is similar to the Hubble expansion rate at the end of the tunneling, we can assume 
the thermal equilibrium. This requirement translates into a condition on hg which is 
easy to satisfy. 5 

Since the warp factor of the S-throat is much smaller than that of the A-throat, 
the KK spectrum is much denser. After the thermalization, the KK modes will decay 
to the much lower mass levels and become much more difficult to tunnel back out. 
Therefore the tunneling from A to S is effectively a one-way process. This streaming of 
the KK modes heats the S-throat. At the end of the tunneling process, the universe will 
become dominated by the thermalized KK modes in the S-throat and become radiation- 
dominated again. 

r ^ + 3 -, + 2 

5 Ignoring all other factors, this condition requires the warp factor hs < r t 5 J^ 16 • mass 5 t+ 16 . For 



example, according to the discussion in Sec. 8.2 and the result flS.q), this is naturally satisfied 
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But there is a subtlety here. Due to the same reason, i.e. the S-throat has a much 
smaller warp factor, the initial thermalized gas may become stringy. To get an estimate 
on the temperature using the KK spectrum (|6.2|), after tunneling we have 



tfun (hsR S ) 7 

Borrowing the result (|8.6| ) from Sec. §, we see that 



2 , r>-L.. • (6-15) 



7+2 1 2 7+2 t+2 



Tsit tun ) < V S +A N£ +A M£ A R S ^ A oc h} +A m s . (6.16) 

On the one hand, we see that this upper bound is typically much smaller than the mass 
level h^R^ 1 of the lowest KK mode in the A-throat; on the other hand, we see that it 
is possible for Ts{t tun ) to exceed the red-shifted string scale h s m s in the S-throat. To 
proceed, we will ignore the involvement of the possible string states for the following 
reasons. First, as universe cools down the thermalized gas will eventually become a gas 
of KK particles. So the estimates on the KK relics still apply. Second, in the most 
pedagogical sense, since ( |S.16|) is just a upper bound, we may adjust various parameters 
to make T s (t tun ) stay below the red-shifted string scale. However, a better understanding 
of such a possible stringy phase and its consequences is clearly very interesting, because 
it may not happen too early away from the BBN. 

Before continuing our story in the S-throat, we note that there is another possible 
period during the tunneling that deserves attention. As we have seen from Sec. |2], if 
bulk resonance happens, KK will stay in the bulk in this intermediate period, so it is 
also important to analyze the graviton production and possible KK relics in this case. 
It is easy to check that any KK in this resonance period is decoupled from each other. 
So any stable KK, once produced through decay process, will be a problem. But this is 
not a big concern, because there is usually no good candidates for stable KK modes in 
the bulk. The wavefunction of these KK modes spread over the bulk, and will always 
hit some throats. This is in fact how we estimated the tunneling rate to X-throat in 



Sec. |23| For graviton production, we integrate all particles and get a very small ratio 

n grav 

P — « g 1/4 Ml[ 2 a Kg T]H , (6.17) 



where T tun ~ h A R A is the tunneling rate from the A-throat to the bulk. 
6.3 KK evolution in S-throat 

The evolution of KK modes in S-throat is similar to that in A-throat. A small difference 
is the additional relativistic degrees of freedom g due to the Standard Model particles 
at the end of the mass spectrum. So we take g to be n y KK if g > go, and approximately 
constant g Q later on. 
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The relative abundance of the possible KK relics is similarly obtained, 



/ n KK\ 

« £ ^^Ulr • (6-18) 

\ P J RDMD riKK 



Comparing to (|6.7|) this is typically negligible since okk here is enhanced by powers of 
h s < h A . 

The graviton is mainly produced around t tun , 

~ 4 /4 (hsTl'f* . (6.19) 

\ P / BBN or RDMD V 7 

For small 7, this ratio is much smaller than ( |6.9|) due to the dilution of the KK particles 
during the tunneling. 

Overall, we summarize the graviton and KK abundance obtained in this section and 
Sec. |] for the double throat case. 

For the graviton, the total abundance is a summation of those from the initial decay 

in Sec [|, (-f) [^~) (where the last extra factor is added due to the matter- 
dominated phase during the tunneling), the annihilation in the A-throat Eq. (|6.9|) , the 
S-throat Eq. (|6.19 ), and the bulk Eq. (|6.17 ) if resonance happens. All of them are 



suppressed by various combinations of warp factors, bulk size and matter-dominated 
phase, so easily become very small. The actual value is parameter dependent. 

For the KK relics, Eq. ( |6.18|) is suppressed by powers of the S-throat warp factor, 
so the A-throat contribution ( |6.7| ), if stable, is typically larger. Evaluating ( |6.7| ) using 
( |3TT9D and (|6l|) , we get 



SIIdmd ~ ( h A M 3 p { 2 ] p- R T MD R A T}l 2 n . (6.20) 

For example, if M P1 w 10 18 GeV, p% MD w 1 eV, m s w 10 10 GeV and N A w 10 4 , for the 
resonance case we need 



R A 

T 



h 5 / < IO- 31 (6.21) 



to satisfy ^^md ~ 1- This for example can be satisfied if L w 10 3 R A and h A ~ 10 -3 . 
Note that the resonance condition typically requires a large L. For the non-resonance 
case, we need 



*± ) < 10- 31 . (6.22) 



This for example can be satisfied if L w 10 R A and h A ~ 10 . 
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7. The bulk case 



Before going to the more general multi-throat case, let us look at the case where either 
SM branes or brane annihilation are in the bulk. 

So far we have considered SM branes in throats. These throats can be regarded more 
generally as black hole like geometry, created by large number of D3-branes, 3-form fluxes 
or even induced large D3-charges on the D7-branes. To enter these geometries KK modes 
need to tunnel. But it is also possible that some small number of branes appears in the 
bulk rather than hiding in some throats. For these branes, the warping they induce 
is small so KK modes do not have to tunnel a significant barrier. In this situation, if 
the KK wavefunction appears in the bulk, it can intersect with these branes and the 
KK particles can decay into light brane fields. So they start to heat the bulk branes 
instead of heating a throat through tunneling. This can happen in the bulk resonance 
case when KK propagate in the bulk, or in the non-resonance case where KK also has a 
wavefunction in the bulk although it is not propagating. Namely these small number of 
branes in the bulk may become SM branes. 

The heating in this section can be easily studied by collecting a few results in Sec. |^ 
and |6]. There are three possibilities. 

We first discuss the case in which the brane annihilation still happens in a throat 
(A-throat). The heating of the SM branes is caused by the KK particles tunneled out 



of the A-throat. The analyses in the A-throat are the same as that in Sec.SJ.. The 
difference is that, the KK particles will tunnel to the bulk and decay to lower KK modes 
and SM brane fields. Same as we discussed at the end of Sec.|6.2|, any stable KK will 
be dangerous in this case, but the requirement of no stable KK is reasonable for the 
bulk since they will always intersect with the bulk SM branes and decay. The graviton 



production is small and the same as Q6.17 ). So the thermal history is rather similar to 



the double throat resonance case and the heating works. 

The second case is when both brane annihilation and SM branes are in the bulk. 
The analyses for the KK modes are the same as in the single throat case. We get the 
same expressions as in Eqs. ( |5.10|) and ( [5.14 ) . But in this case, there is no longer warp 



factor and bulk size enhancement factors to the KK interactions, &kk ~ M p f and we 
get 

^RDMD ~ Yl 9~ 1/4M PlPRDMD • (7-1) 
riRK 

However we still do not worry about the KK relics in the bulk, because they are unstable 
against decaying into brane fields. The graviton production during the KK evolution is 
still very low. Overall, for this case, the main graviton abundance problem comes from 
the initial brane and string decay period discussed in Sec.^. Namely, the analyses in 
Ref. shows that in flat space the graviton is produced more than KK particles in the 



initial string decay. In all other cases that we considered, the brane annihilation happens 
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in a throat so this graviton problem is absent. This is the only case that can have a 
serious problem in heating. This problem is present in models such the branes-at-angle 
56], the D3-D7 inflation with D7-branes in the bulk |57|, or the D3-anti-D3- 



sccnario 



brane scenario in the bulk |2|, |3|] , if the SM branes are also in the bulk. 

This graviton problem can be cured if there is a matter-dominated era before the 
BBN, so the graviton density can be considerably red-shifted. This leads to the third case 
where we locate the SM branes in an S-throat, although the brane annihilation happens 
in the bulk. If there is no extra small number of (anti)branes left in the bulk, the KK 
modes will eventually tunnel into the S-throat. This introduces a matter dominated 
phase similar to what we have seen in the double throat case. The graviton problem is 
hence cured by the red-shift and the story of the S-throat in Sec. |6.3| applies here. 

8. Multi-throat heating 

In this section, we consider the heating in a more general multi-throat compactification. 
In this case, the brane-anti-brane annihilation at the end of the inflation still happens in 
the A-throat. The heating from the A-throat to the bulk or the S-throat follows the same 
analyses in the previous section, but the presence of many other throats gives additional 
interesting phenomena. 

8.1 Throat deformation in Hubble expansion 

The spectrum in a throat is mainly determined by its minimum warp factor. However in 
the Hubble expansion background, the warping of a throat cannot be arbitrarily small. 
Even if we start with an infinitely long throat, the Hubble energy will impose an infrared 
cutoff on the warp factor. 

This Hubble back-reaction may be analyzed in different ways |58|, [37], [59[ . One 



way to look at it is to consider the back-reaction of the closed string created in the IR 



side of the warp space ||12|| . Assuming an infinitely long throat, such string creation 
will be inevitable because the local string scale will be red-shifted down the throat and 
eventually fall below the Hubble energy scale. The energy density of these strings is H 4 . 
The back-reaction of such closed strings become significant to the source of the throat if 

H 4 » hi m NT 3 , (8.1) 

since the throat with R A « N/T 3 can be thought of being sourced by N number of 
branes. Therefore, although in static background we can have multiple throats with 
arbitrary warp factors, in the Hubble expansion background they will have a cutoff at 

HR 

v N 

if they are too long. 
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Another way to look at it is to consider various moduli that stabilize the throat |37| . 



Although in static background, throat with arbitrary warp factors can be constructed by 
stabilizing those moduli, the moduli mass will obtain a correction in a Hubble expansion 
background, and therefore cause the deformation of the throat. The explicit evaluation 
along this line has so far been difficult due to the subtleties analyzing the supergravity 
equations in the presence of small warp factors Q . 

Here we give a toy model of this type to roughly estimate this effect. Consider an 
infinitely long throat generated by N stack of branes rather than fluxes. In the Hubble 
expansion background, the brane will fluctuate in the transverse directions, 

Acf) « v^Ar « ^ . (8.3) 

This fluctuation causes a deviation of the brane moduli which otherwise are located 
exactly at the origin. So the throat is no longer extended to infinity at r = but 
terminated at a minimum r min « Ar. This corresponds to a minimum warp factor 



We see that in this toy model we get a consistent picture as in Eq. (|8.2|) . This consistency 
can be regarded as a duality relation between deformed AdS and CFT: Eq. ( |8.2| ) is 
obtained from the point of view of the gravity side, while Eq. fl8.4j) the field theory side. 

These analyses apply to both the inflationary and heating epoch, and give a time- 
dependent minimum value for throat warp factor at corresponding moment. As a remark, 
we note that the energy associated with those closed strings and moduli potentials is 
negligible to the heating energy, since it is of order 0(H A ) and much less than the 
inflationary energy. But it can create a spectrum of cosmic strings with relatively low 
tension Ill2fl. 



8.2 Selection of long throats 

In the presence of many throats, the tunneling of KK modes can have many channels. 
That is, it is possible for the KK modes to go from A-throat to different throats. So 
which throats become dominantly heated depends on the tunneling branching ratio. 
A necessary condition for the tunneling discussed in Sec. |2.2| and |2.3| to happen is 



the matching of energy levels in both throats. This condition requires that an energy 
eigenvalue in the X-throat fall within the energy width of the initial state in the A-throat. 
The width of the lowest KK mode in the A-throat is determined by its tunneling rate 
T tun , as one can see from the uncertainty principle. So, the more finely spaced is the 
X-throat spectrum, the more likely such a matching can be found. 6 A long throat is thus 
6 Note that this is different from the argument that the large phase volume in a long throat guarantees 
more heating after string decay because of equipartition, as discussed in Ref. ^5| . As pointed out 
in the same references, the presence of the potential barrier invalidates such arguments. Here even for a 
short throat, if a mass level matching is accidentally satisfied, the tunneling will still proceed (but may 
oscillate). Decreasing warp factor increases this chance. 
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favored during the tunneling heating. Since in a throat with a minimum warp factor hx, 
the energy spacing for the s-wave is hxR-x > if 

h x < T tun R x , (8.5) 

the tunneling will generically happen for this throat. The Hubble constant during the 
tunneling is H w « so the allowed minimum warp factor ( |8.2| ) in this Hubble 
expansion background is 

h min « . (8.6) 

V-'Yx 

So the condition (|8.5|) is allowed. 

Thus we have seen a dynamical mechanism of selecting a long Standard Model 
throat in the multi-throat heating. For example, for the tunneling rate ( |2.32| ), if roughly 
Ra ~ Rx ~ D, Ttun ~ h A D~*, any throat with 

hx < h\ (8.7) 

will be generically heated, while heating probability for a shorter throat is smaller. A 
throat with a large hierarchy ratio (~ 10~ 10 ) in the sense of Randall and Sundrum can 
be naturally selected in this way. If the SM branes are in the bulk, such a mechanism 
will favor a large bulk size. 

As a comment, we note that, even if the mass levels between the A and X does not 
match, tunneling can proceed through a secondary effect. A virtual KK particle is first 
created in the X-throat, and then splits into several lower mass states such as lower KK 
modes or brane fields. The kinetic energy of these light particles can compensate the 
mass level difference. It will be interesting to estimate the branching ratio of such a 
process and we leave it for future study. 

8.3 Warped KK dark matter and hidden dark matter 

In this subsection we discuss two interesting possibilities of non-relativistic dark matter. 
The first is the KK relics in the S-throat. In our setup, the lowest KK mode in the 
S-throat will interact with the SM branes and decay into brane fields, i.e., SM particles 
as light open string modes. However, higher level KK modes can have a long enough 
life time, e.g., if they are associated with a conserved angular momentum in the throat. 
For example, in the Klebanov-Strassler throat, besides the excitation quantum number 
n, the KK modes can also have angular quantum numbers j for the approximate S 3 and 
/ for S 2 in a realistic throat, i.e., a warped deformed conifold. The (approximate) S 3 
remains finite at the bottom of the throat (a property of the deformation of the warped 
conifold). As a consequence, j > 1 quantum number costs less energy than the I quantum 
numbers ||40|| . The KK modes with only j quantum numbers have non- vanishing wave 
functions at the bottom of the throat, so they can easily decay to SM particles. (Recall 
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that transverse momenta can be absorbed by the solitonic-like branes.) On the other 
hand, the S 2 shrinks to zero at the bottom of the throat, so the wave functions of I > 
KK modes are suppressed at the origin, much like that in the hydrogen atom. Having no 
overlap with the SM branes, the light / > modes will not decay to SM particles. They 
couple to SM particles via the exchange of a s-wave KK mode, so two I > KK modes 
can annihilate. However, this annihilation rate drops rapidly as the universe expands, 
and these KK relics abundance can be estimated as we did for the simplified throat in 
the previous sections. Their tunneling rate is very small due to the smallness of the 
S-throat warp factor and additional non-s-wave suppression factors. So these lightest 
angular KK modes are quite stable. 

Such a non-relativistic relic was a concern in Ref. |35] since their abundance was 
estimated to be extremely large. So if stable they would over-close the universe. In the 
previous sections, we have analyzed the KK thermal evolution and the abundance of 
these KK relics more carefully. As we discussed, for example in Eq. (|5.5|) , these relics 
typically decouple non-relativistically. From Eq. ( |5.11| ) and ( |6.18|) , we can see that the 
relative density of these KK relics is suppressed by warp factors, bulk size and matter- 
dominated tunneling phase, and so can be naturally made safe. (To compare with the 
analysis in Ref. |ft5| , we note Eq. (4.16) in JJIJ is based on relativistic decoupling and 
has no adjustable small parameters. Here Eq. ( p. lip and Q6.18] ) have several important 
differences: a factor of h 2 AS due to the non-relativistic decoupling; a factor of (R/L) 6 
from the throat- versus-bulk effect; and in the double throat case, an extra suppression 
(tAnr/Uun) 1 ^ 2 for the A-throat from the matter-dominated tunneling phase.) 

Of course, one can then tune various parameters in the formulae to make them just 
account for the density of the observed dark matter. Namely they can become the warped 
KK dark matter. In the single throat case, the relics abundance is given by ( |5.11| ); in the 
double throat case, ( |6.18|) is typically smaller than ( |6.7|) , but it may become important 
if the KK relics in A-throat have too short (tunneling or decay) lifetime. (Otherwise we 
will discuss it in Sec.|9.2|.) They may be detected by collider experiments or underground 
dark matter searches. For example, when two SM particles collide at the red-shifted 
KK energy scale, they can produce an intermediate unstable s-wave KK particles, then 
subsequently decay into two warped KK dark matter particles. They will have different 
cross-sections than those of the more standard dark matter candidates. Discussions on 



various other kinds of KK dark matter can be found in e.g., Ref. [41, 42, 43, 44, 45 



Next we consider another type of dark matter. In the multi-throat setup, when the 
KK particles are tunneling from the A-throat to heat the S-throat or the bulk SM branes, 
similar process can also happen for any throat which satisfies the mass level matching 
condition discussed in Sec. |0| and Sec. Q In particular, such a tunneling happens 
generically for any throat long enough to satisfy ( |8.5p . Even in the single throat case, 
it is possible for a tiny amount of KK modes to tunnel out to other throats. 7 In this 



7 The possibility that KK particles may tunnel to other throats is also discussed in Ref. [B5l p6[ 
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subsection and Sec. PTj we propose that the energy density going into those throats 
can create new dark matter candidates. Using the knowledge of the thermal evolution 
that we obtained in this paper, we will discuss the properties and thermal histories of 
these dark matter candidates and the plausibility of this proposal, in particular paying 
attention to the dark matter coincidence problem. We call those throats as D-throats, 
with D stands for dark matter, and call those dark matter as hidden dark matter, because 
they are in a hidden sector and their wavefunction has exponentially small overlap with 
our visible sector. 

The identities of the hidden dark matter in D-throats can be very different from 
those of the warped KK dark matter. For a D-throat without any branes in the end, 
as in Sec. [Q] , the KK tower will eventually be lowered to the ground state. So the 
hidden dark matter in such a throat will be the lowest KK particles, but not those 
with conserved angular momentum. These KK particles are mutually decoupled and 
cannot find each other to annihilate. Different from Sec. |6.1| , since the D-throat has a 
much smaller warp factor than the A-throat, these lowest KK modes take much longer 
to tunnel out and therefore become stable. For a D-throat with branes in the end, the 
identity of the hidden dark matter can be either similar to that of the warped KK dark 
matter in the S-throat, or more exotic particles analogous to our baryons. 

These dark matter have very different properties comparing to the usual dark matter 
candidates such as the warped KK dark matter mentioned above, the SM KK dark matter 



|44|, [45| , axions |K| ^T], or the lightest supersymmetric particles (LSP) |]63], |65 



1. Coupling: For the usual dark matter candidate, although they are decoupled to 
the Standard Model particles at low energy, they can couple if the energy is excited 
reasonably higher. For example, for the warped KK dark matter in the S-throat, 
the corresponding energy scale is the warped KK scale, lower than the red-shifted 
string scale; for the LSP in the MSSM, it is around the supersymmetry breaking 
scale TeV. 

However, for the hidden dark matter in the D-throats, it almost completely decou- 
ples from the Standard Model. As we can see from Eq. ( |2.17|) - (|2.19|) , the wave- 
function of a KK mode damps at least by a factor of h 8 D from D to S-throat. The 
particles in the S-throat can couple to the hidden dark matter only if they are ex- 
cited out of the throat, namely by roughly the UV string scale. Hence practically, 
the only interaction between the hidden dark matter and the Standard Model (or 
its extension) is through graviton mediation. Of course the hidden dark matter 
can have self-interactions. 



Ref. |3J| raised a concern that this would generate KK relics that over-close the universe and destroy 



the BBN. In Sec. 9.1, we will see how the unusual properties of the hidden dark matter can be used to 
avoid such problems. Ref. |?6| gave a lower bound on these throats to avoid the overheat of the S-throat 
after BBN due to KK tunneling from those throats. 



-37- 



2. Thermal history: For the usual dark matter, they initially thermalize with the rest 
of the particles and closely participate in the thermal history when the temperature 
in the S-throat is high. Only later on after freeze-out, they evolve independently. 

However, the hidden dark matter in D-throats has its own thermal history. They 
have different temperatures once particles thermalize in the D-throats. The thermal 
evolutions of the particles in the D-throats and the S-throat are independent, with 
the net effect driving the spatial expansion of the whole universe. 

3. Relic density: For the usual dark matter, the relic density is set by the decoupling 
condition in the course of the thermal evolution. The decoupling happens when the 
density of the particles drops low enough so that their annihilation rate is below 
the expansion rate. After that, they are driven apart forever by spatial expansion 
and cannot find each other to annihilate. We have seen such an example in the 
discussions of the KK relics in the previous sections. 

But for the hidden dark matter in D-throats, it never evolves in the thermal history 
of the S-throat and decouples at the beginning as soon as KK particles tunnel to 
the D-throats and decay. The relics density is set by the initial conditions during 
the tunneling and the properties of the D-throats. 



9. Hidden dark matter 



In the first subsection, we study more details of the hidden dark matter in D-throats 



proposed in Sec. |8]^ and the dark matter coincidence problem. In the second subsec- 
tion, we discuss another type of hidden dark matter in the A-throat and its possible 
implications on cosmic rays. 

9.1 Hidden dark matter and dark matter coincidence problem 

We first list some interesting consequences and constraints followed from the properties 
discussed in Sec. E73J In this subsection, when we make explicit comparisons, we will 



choose to compare the hidden dark matter with the LSP dark matter for example. Also 
we will choose S-throat as the location of SM branes. The discussion still applies if we 
replace it with SM branes in the bulk. 



1. The mass of the LSP has to be greater than 100 GeV, because otherwise it would 
have already shown up in colliders. But this is not the concern for the hidden dark 
matter. In fact there is little constraint on the mass of the hidden dark matter 
particles. A more relevant statement is that the hidden dark matter does not have 
to become non-relativistic before the temperature is 100 GeV in the S-throat. It 
only has to take effect before the radiation-matter transition epoch in the S-throat 
at around 1 eV. 
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2. Each throat, including the S-throat and D-throats, has its own RDMD transition 
epoch. To be consistent with our observation, the S-throat has the latest RDMD 
transition. The RDMD transition in D-throats can have different nature, depend- 
ing on their matter contents. We will give two examples later. For an observer in 
a D-throat, after its local RDMD transition, she will observe a period of hot dark 
matter-dominated era, contributed by the S-throat. 

3. From the observations, we know that the density of the dark matter is approx- 
imately 5 times of that of the baryons. If the hidden dark matter in D-throats 
accounts for most of these dark matter, we require 



n Z(+RDMD 

n { .RDMD\ a \ l D )_ ^ k f.RDMD\ / Q -, \ 



RDMD\ a \ Z D ) k /.RDMD 



where the subscript D or S refers to the local quantity of the D or S-throat respec- 
tively. 

Moreover, at the time of BBN, we require the density (of either matter or radiation) 
in the D-throats be much smaller than that in the S-throat, so that the process of 
the BBN is not affected. Therefore during the tunneling heating, the most heated 
throat will turn out to be the S-throat. This is consistent with the previously 
mentioned requirement that the RDMD transition in the S-throat happens last. 
In fact, the energy distribution between the S and D-throats does not have to be 
too different. The requirement 



PD(t 



BBN\ 
S > 



Ps(t 



BBN^ 



< a few percent (9.2) 



will be enough to be consistent with the observations. This variation of the ini- 
tial densities is due to the variation of the tunneling rate with respect to different 
throats. Indeed from the formulae (|2.31|) , (|2.32|) and (|8.5|) , we can see that, while 



the tunneling rate has the same scaling dependence on the warp factor tiA, varia- 
tions are expected for D-throats having different locations, length scales or warp 
factors. 



We now look at two examples with different contents of hidden dark matter in D- 
throats. We will focus on the following coincidence problem: why the density of the 
dark matter happens to be roughly the same as (five times) that of the baryons? For the 
usual dark matter, this density is formed when the density of the dark matter particles 
exponentially decreases and eventually freezes out after the annihilation rate drops below 
the Hubble expansion rate. Then this question amounts to showing whether the required 
annihilation cross section is naturally obtained given a model. For example, it has been 
argued that the WIMP may naturally have this cross section. Here, since we have very 
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different properties, it will be very interesting to see what new aspects hidden dark 
matter offers to this classic question. 

• Example 1: Imagine a D-throat without any brane located at the tip, so the hidden 
dark matter will be the lowest KK mode with mass m^ K ~ hpFC^. This is the case of 
hidden KK matter. 

In this example, the RDMD transition in the D-throat happens when the local 
temperature crosses m^ K at Tjj(t^ DMD ) ~ hpR^ 1 . At this moment, the coincidence 
problem requires the energy density in the D-throat to be five times of the baryons in 
the S-throat. Since the total energy density is still dominated by radiation in the S-throat 
until tg DMD , we obtain the ratio of the energy density between the D and S-throat at 



and before t% DMD ', 



„ T (+RDMD\ 

PR „ } - (9 3) 

n T (+RDMD\ ■ y J -°> 

Ps J-sytn ) 



Suppose initially we have (|9.2|), Eq. ( |9.3| ) indicates that 

t RDMD > 1Q 4 t RDMD (g ^ 

To get a more explicit requirement on the D-throat, we use Eq. (|6.2|) and find initially 

(9.5, 

T s \PsJ \h s R s L J 
At the local RDMD transition in the D-throat, 

T D (t* DMD ) « h D R~ D l . (9.6) 

At the same time, we know from the condition ( |9.3| ) and Ts(t§ DMD ) ~ leV that 



T s {t* DMD ) « ^ • 5 eV . (9.7) 



Pd 



Plug Eq. (pD and (Q into QOD , we get 



7+3 



Pd J \R 



Ps\ 4 (Rs\ i (5 eV RD ^ h -l (gi 



In this example, the coincidence problem results in a tuning of the D-throat warp factor 
for any ratio ps/pd, as long as the condition (|9.4|) is satisfied. For the bulk SM branes, 
one simply replaces ps with pbuik, hg Rs with L in Eq. ( |9.8D . 

• Example 2: Perhaps the S-throat does not have to be special in the sense that 
branes can also be located in other throats such as D-throats as well. In each D-throat, 
there are different versions of the "Standard Model" on branes. In such a case, the 
relevant hidden dark matter content is no longer the lowest KK mode, it is different 
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matter fields on branes, analogous to our baryons. More importantly, the mass of these 
matter is no longer of the local KK mass scale. The hidden dark matter is massless at 
the leading order comparing to this scale, just like the baryons in our Standard Model 
throat. This is the case of hidden matter. 

Again we naturally assume that S-throat and D-throats have similar but variable 
heating branching ratios (for example, if they all satisfy (|8.5| )). We pick the most heated 
throat as the S-throat for reasons discussed above. The thermal evolution of the D- 
throats are similar to the S-throat. As the universe expands, the KK tower lowers and 
eventually transfer its energy to the brane fields. The branes fields are mostly massless 
particles, but there can be a tiny fraction of asymmetric "baryons" matter generated 
by "baryogenesis" . Since radiation red-shifts faster than matter, the "baryon" matter 
eventually dominates in the D-throats. What follows is a period of "hot dark matter" 
dominated era contributed by the S-throat. After that the whole universe becomes 
matter-dominated. 

In this case, the coincidence problem reduces to the similarity between the "baryoge- 
nesis" efficiencies on the "Standard Models" in the D-throats and that on the Standard 
Model in the S-throat. To give an example, we assume initially ps ~ 20p/> If the D- 
throat matter density is similar to ours so that p™ atter ~ 0.2p^ atter , the RDMD transition 
happens a little earlier in the D-throat - by a factor of 100 in terms of the scale factor. 
That is, Eq. ( |9.4j ) still applies here. By the time that the RDMD transition happens 
in the S-throat, we have both the baryons and hidden matter with similar abundance. 
Both the mass of the hidden dark matter and its final abundance become much more 
independent of the D-throat warp factor. 

Both examples give very interesting illustrations of how hidden dark matter can 
behave differently against the usual intuition. Naively, if the energy density that goes 
into a dark matter candidate (e.g. LSP) is as large as a few percent of that goes into the 
SM particles, it surely over-closes the universe and destroys the BBN. This is because the 
LSP dark matter candidates become non-relativistic before the temperature is 100 GeV. 
In terms of time, this happens more than 10 22 times earlier than the RDMD transition. 
Only extremely tiny fraction of energy is allowed to deposit into the dark matter to 
account for such a large red-shift. But hidden dark matter can become non-relativistic 
much later, as late as in (|9.4j ) - reducing 10 22 to 10 4 - which greatly shortens the red- 
shift. This is achieved, in the first example, by having a long D-throat and therefore very 
low KK levels; in the second example, by having similar amount of matter and radiation 
between the hidden sector and the visible sector. 

9.2 Hidden dark matter and cosmic rays 

So far in this section, we have studied the properties of the hidden dark matter lying 
in D-throats which are different from both S and A-throat. But in the double throat 
case or the bulk case, the KK modes with specific conserved angular momentum in the 
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A-throat can also become viable dark matter candidates. It is possible that these KK 
modes have very small tunneling rate. For example, the s-wave in the A-throat have a 
tunneling rate h\R~^ or h^R^ 1 - This s-wave is responsible for the tunneling heating. 
For a l-th partial wave, the tunneling rate is expected to be reduced by a factor of h\ 
for the simplified geometry studied in Sec. ^| (same as the absorption rate in [53[]). Hence 



a reasonable higher partial wave can be stable against the tunneling. If this angular 
momentum is also approximately conserved (now it can be associated with either S 3 or 
S 2 in the KS throat since we do not have branes left in the A-throat), it becomes a 
candidate of stable relics. We have estimated the relics abundance of such KK modes in 
Eq. (|6.7| ). We see that it again contains various tunable small parameters and can be a 
viable dark matter candidate. 

An interesting feature of this dark matter is that its properties lie in between the 
hidden dark matter in D-throats and the warped KK dark matter in S-throat. In terms 
of couplings, it is similar to the former — it only couples to the visible sector through 
gravity and usually cannot be produced in colliders. So it is still the hidden dark matter. 
But in terms of the properties of the thermal history, it is more like the latter — it 
decouples as a non-relativistic thermal relics and its abundance can be calculated in the 
usual way. 

We have been familiar with the fact that large cosmic strings can be left in the 
A-throat as a result of closed string evolution. For these cosmic strings, they cannot 
tunnel. The evolution of cosmic strings will lead to a scale-invariant spectrum (i.e., 
cosmic strings will come in all sizes in some distribution). In general, that will yield a 
density depending on string tension /x. It goes like G^iiT of the critical density, where 
T is a numerical number that is around 10 to 100 f64|. With Gjyfi < 10~ 6 , large cosmic 
strings contributes very little to the critical density. But now we see that it is also 
possible that the closed string and KK evolution lead to some stable KK modes in the 
A-throat. Then those KK modes can have a sizable density to contribute to the dark 
matter. 

Such a dark matter is typically very heavy. Take the double throat non-resonance 
case as an example. In the numerical example of Eq. (|6.22|) , the mass of this hidden 
dark matter is 

M DM ~ 10 7 GeV . (9.9) 
One can easily change the string scale since the bound ^r^ md < 1 on fl6.20| ) is not 



very sensitive to this variation. For example, taking m s « 10 13 GeV, Na « 10 2 and 
still keeping the tunneling rate T tun ~ h l A 7 R~]^ and the other parameters same, we get 
Ha ~ 10 -2 and L ~ 60R to saturate the bound. In this example, the mass of the hidden 
dark matter becomes 

M DM ~ 10 n GeV . (9.10) 
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There are many interesting consequences of such a dark matter. 

First, we have seen that this hidden dark matter belongs to thermal relics, so both 
( |9.9|) and ( 9.1 0|) evade a well-known bound |)C] which states that the maximum mass of 



a dark matter is a few hundred TeV if it is a thermal relic. This bound comes because 
larger dark matter mass leads to smaller annihilation cross section and more red-shift, 
and therefore excessive dark matter density. The reason that we evade it is the extra 



matter dominated tunneling phase discussed in Sec. |6.1| and |6.2| . During this phase the 
total energy of the universe, including that of this dark matter, is hung up in the A- 
throat in a long matter-dominated era. The radiation-dominated phase in the S-throat 
is therefore much delayed. This greatly increases the allowed ratio between the density 
of the dark matter and the rest. 

Second, this dark matter can have a lifetime longer than the age of the universe, 
but occasionally tunnel to the visible S-throat. Namely this becomes a novel example 
of an ultra heavy thermal relic being the source of the ultra high energy cosmic rays. 
To estimate the lifetime of those angular KK modes in the A-throat, we use the ansatz 



h A Ra- In the above example ( |9.10D , we can see that the lifetime is longer than 



the age of the universe 10 17 s if Z > 2. 

Third, since this cosmic ray is due to the decay of the dark matter in our galaxy, they 
do not have to travel a long distance between galaxies to lose energy through interactions 
with CMB. So the GZK cutoff |67], ^] does not apply. The energy of these ultra high 
energy cosmic rays can be as large as 10 20 eV such as that in Eq. ( |9.10|) . Since the source 



of such decay is not local objects, the distribution of the cosmic rays is more isotropic. 
Comparing to the distribution of the dark matter annihilation, which is proportional to 
p|, M (x), these cosmic rays will have a distribution proportional to /3dm( x )- Moreover, 
the energy of the cosmic rays Mqm correlates with the inflation scale Vmf and the tension 
of the cosmic strings /i^, 

M DM « n KK {g s N A )- x I^T « n KK N- 1/4 V^ . (9.11) 

If N^ 4 is not too big, they are all close to each other in terms of order of magnitude. 
So in this scenario a measurement of the energy of the relevant cosmic rays predicts the 
primordial gravitational wave and the cosmic string tension. 

Fourth, the nature of this decay can be extremely stringy. This happens in the 
double throat case. Although the original KK mode in the A-throat is a field theory 
mode, after it tunnels to the S-throat, its mass is way above the local red-shifted string 
scale due to the large warping of the S-throat. So detailed observations of such decays 
may become a useful laboratory to compare with the theoretical calculations in string 
theory. 

9.3 Comparison to other dark matters 

Besides the LSP that we have discussed, in this subsection, we give a couple of more 
comparisons between other dark matter candidates and the two main categories of dark 



- 43 - 



matter candidates proposed here, i.e. the warped KK dark matter and the hidden dark 
matter. 

The warped KK dark matter shares common features of most usual dark matter 
candidates, such as their thermal properties and coupling to SM particles. For example, 
conceptually it is similar to that proposed in the 5-d RS geometry [(41], [4^]. Besides the 
differences caused by those between the 5-d RS and 10-d GKP setup, there is no graviton 
KK in 5-d stable against decaying on the IR branes, so in jll, |2] , such a stable particle 
was put in by hand ("bulky"); here we have a natural and viable realization of their 
"bulky". 

The hidden dark matter is more unusual. Our definition of the hidden dark matter 
is that they are in a hidden sector, and their wave-function only have an exponentially 
small overlap with our visible sector. As far as detailed properties are concerned, they 
may have overlaps with various other kinds of dark matter candidates. But they also 
have differences. For example, superWIMPs or GIMPs (super-weakly or gravitationally 
interacting massive particles) such as gravitino or graviton KK modes in universal extra 
dimension (UED) also have gravitationally suppressed interaction with other particles 
|S9f1 . But these particles, for example the gravitino, can be the decay product of a su- 



persymmetric Standard Model particles with a mass such as 1 TeV. Hence superWIMPs 
actually couple stronger than the hidden dark matter (but weaker than the warped KK 
dark matter in the S-throat). For example, if one collides two SM particles at 1 TeV, 
they can become an intermediate supersymmetric particles, through a relatively strong 
coupling, and then decay to gravitino, through gravitationally weak coupling. But for 
hidden dark matter, even if the energy of the colliding particles is of order their mass 
level, they still can only be created through the mediation of gravitons. Because of the 
same reason, the gravitino can inherit the usual relics thermal abundance of the mother 
particle, but hidden dark matter in the D-throats has a completely independent thermal 
history. As another example, there is a proposal of having ultra high energy cosmic rays 
as decay products of non-thermal heavy dark matter ]70| . But as we have discussed, our 
heavy hidden dark matter in the A-throat is thermal relics, while the non-thermal hidden 
dark matter in the D-throats does not have to be heavy at all. Of course, more detailed 
study and comparison with various dark matter candidates is necessary and interesting. 



10. Conclusions and discussions 

In this paper, we have studied the properties of the graviton KK modes in the single 
and multi-throat compactification, including their wavefunctions, tunneling properties 
and interactions. We have studied the thermal history after the brane annihilation 
at the end of the brane inflation, to test whether the heating is compatible with the 
standard big bang nucleosynthesis. We have found that, as long as the brane annihilation 
and the SM branes are not both in the bulk, the heating process is viable no matter 
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where they are: one of them in the bulk, in the same throat or in different throats, 
although different scenarios have different thermal histories. The KK thermal relics and 
the graviton abundance are suppressed by various small factors: warp factors, throat 
versus bulk size, and the red-shifting factor during the tunneling period. The long 
tunneling phase required in some cases is not only compatible with a viable heating, but 
also gives rise to several interesting phenomena: an extra suppression for the KK relics, 
a dynamical mechanism of selecting a long Standard Model throat, and the generation 
of hidden dark matter. Our results on the KK thermal history and relics abundance 
show that there are three types of new dark matter candidates: the warped KK particles 
with specific angular momentum in the S-throat; the hidden dark matter in D-throats; 
and the hidden KK particles with higher partial wave in the A-throat. Physically they 
differ from each other in terms of either the relics thermal history or couplings to the 
visible sector. The novel properties of the hidden dark matter have especially interesting 
implications on the dark matter coincidence problem and ultra high energy cosmic rays. 

We have seen several possibilities that stringy excitations are produced: in the A- 
throat during the initial brane decay, in the S-throat after the tunneling from A to S, in 
any throat with too small warp factors. In the first two cases, the energy density in the 
stringy excitations can dominate the total energy density initially and then lose it to lower 
KK modes. In the last case, it never dominates. It will be interesting and challenging 
to understand if there is any novel effects associated with these phases, besides the 



cosmic strings discussed in Ref. [19], [2(J, ^T], |22fl. For example, it is entirely possible to 



produce some relatively long lived massive open string modes in the SM branes. These 
modes will evolve as a -3 for a while, eventually decaying to the light SM particles. They 
tend to enhance the SM entropy contribution during BBN. However, without an explicit 
model, it is difficult to assess the importance of this possibility. Instead, we see that the 
graviton KK modes play that role beautifully, but in a way that is more intricate. Here, 
the important role played by the warped geometry cannot be over-emphasized. 

The graviton KK modes in different throats are prime candidates as dark matter. 
The stable KK modes in the S-throat have qualitatively similar properties as usual 
dark matter candidates, but with different production and scattering cross-sections. On 
the other hand, the KK modes in the D-throat or A-throat away from the SM branes 
have rather different properties. Since they do not interact with SM particles or their 
superpartners except via gravity, collider experiments and/or underground dark matter 
searches are useless. One way to detect the hidden dark matter is via gravitational lensing 
[7jJ . With better CMBR data, one may fix the parameters in brane inflation more 



precisely. Together with structure formation simulations, rotational curves etc., they 
should allow us to better determine the KK spectrum and interactions, their production 
and subsequent evolution as well as properties of the D or A-throat. One can then check 
to see if their distribution in our universe matches the lensing observations expected in 
the near future. 
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It is quite possible that the hidden dark matter is only a component of the total 
dark matter expected. If direct detection of dark matter candidates and their subsequent 
measurements cannot account for all the dark matter expected, this may be taken as an 
indirect evidence of hidden dark matter. Violation of the GZK bound on cosmic rays 
can be a different piece of evidence. 



The UV DBI inflation [EL 10| and some special case of the IR DBI inflation [11, 12 



require relativistic brane collision as heating, in addition to or as an alternative to the 
brane annihilation. While we expect the main thermal history after thermalization to 
be similar to one of our cases here, it is interesting to see how the details may differ 



especially at the initial moment |f2 



In this paper, we try to give an overall picture of heating towards the end of brane 
inflation. It is an interesting story all by itself. In doing so, we have made only order- 
of-magnitude estimates and neglected all the numerical factors. A more precise analysis 
will probably increase the likeliness of some scenarios and decrease that of others, so 
many interesting aspects are worthy to be studied in more details. It is clear that we 
have considered only some of the simplest scenarios. There are variations of the scenarios 
considered here that deserve investigations. 
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A. a KK cannot decay to all gravitons 

In this appendix, we give a general argument that a KK particle cannot decay to all 
gravitons, generalizing that in Sec. [| 

The 10-dim metric contains a warped four large dimensions and six internal dimen- 
sions, 

d> s io = 9^( x ,y) dx^dx" + g mn (y) dy m dy n 

= A- l ' 2 (y) g^x, y) dx»dx v + g mn (y) dy m dy n . (A.l) 

The warp factor A~ 1 ^ 2 (y) and the internal metric g mn (y) is generated by the matter such 
as fluxes and branes. The wave-functions of KK modes have non-trivial dependence on 
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the internal coordinates, g^ u = r]^ + h^„(x,y), while that of the graviton does not 
depend on the internal coordinates, g^ u = r] pu + h 9 ™ v {x). The 10-d action is 

S W = M( 8 1Q) J d 10 X^—g(XQ)R(XQ) + ^matter • (A-2) 



To get the n-point coupling, we perturbatively expand the action (|A.2|) to the poly- 
nomial of h^ u and d p h piJ . To examine the amplitude of one KK mode coupled to any 
number of gravitons, we replace one of h pv with the wave-function of the KK mode in all 
possible way in each term. We then factor out the wave-function A~ 1 / 2 h^ l f , integrating 
by parts if necessary. This procedure is the same as what we do to get the equation of 
motion, except that we now replace 5g^ u with A~ l l 2 h^ . So we get 

Sio D -Mf 10) J d^X^g^ (A-Wh% K ) - \g» v RW - j^T^A (A.3) 



Under the metric ( |A.l 



< 0) = fljg - \^%)9^ ~ \g pX d m g,ud m g pX + \g px d m g up d m g, x , (A.4) 

where R^} is constructed from 4-d metric g pu , and V(6) is constructed from the internal 
metric g mn . For the background warped geometry, the source does not depend on the 
large dimensions, 

ft /"matter 

:0. (A.5) 



dg 

Therefore, T m / atter is proportional to g ph 



1 



Timatter _ /"matter j \ p\ 

HV — —-^9\w>~' ■ 

Except for the one KK mode, all the other g pv in ( |A.3|) are gravitons and therefore 
y- independent. Hence in the last three terms of (|A.4j ), g pv factors out, and we get 



_ \n P(10) _ 1 ^matter 
Z JW (10) 



g,u [ -^ 2 {6) A- 1/2 - l -A^d m A~^d m A-^ + \a-^R {&) 
+ ^W^A-^C™^ , (A.7) 

/M (10) 



where R$ is constructed from g pu . 
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The last two terms in (|A.7|) cancel because the background equations of motion are 
just these two terms with the overall factor g^ v replaced by r]^ u , the rest does not depend 
on g^ u . The first two terms also cancel because the graviton wave-function satisfies 
R^y — T-gpvR^ = by definition, which is beyond the linear approximation and includes 
all possible non-linear couplings. Therefore ( |A.3| ) vanishes identically, namely, the decay 
product of one KK mode cannot be all gravitons. 
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